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Abstract 


When  coherent  light  is  incident  on  a surface  which 
is  rough  in  comparison  to  optical  wavelengths,  a random 
interference  pattern  called  speckle  is  produced  in  the 
reflected  field.  The  effects  of  speckle  on  an  imaging  system 
composed  of  a laser  scanning  a rough  surface  and  a receiver 
measuring  the  reflected  field  are  examined.  Three  types  of  lasers 
and  two  types  of  receivers  are  considered.  The  lasers  are  designated 
by  their  bandwidths  as  monochromatic,  quasiraonochromatic , and 
broadband.  The  two  receivers  are  the  direct  detection  and 
the  heterodyne. 

The  field  reflected  from  the  rough  surface  is  "crudely" 
modeled  by  multiplying  the  incident  field  by  a reflectance  term 
and  a random  phase  term  which  is  indicative  of  speckle  effects. 

By  second  moment  techniques,  the  imaging  systems  are  modeled 
in  terms  of  their  mean,  covariance  and  system  parameters.  All 
fields  are  propagated  by  the  Huygens-Fresnel  integral. 

The  system  consisting  of  a broadband  laser  and  a direct 
detection  receiver  is  shown  to  be  independent  of  speckle  effects. 

All  other  systems  are  shown  to  have  "noise",  due  to  speckle 
effects,  which  is  signal  dependent.  The  signal  to  noise 
ratio,  for  the  systems  which  include  a direct  detection  receiver, 
is  found  to  be  greater  than  for  systems  which  have  a heterodyne 
receiver.  The  systems  with  broadband  lasers  are  shown  to  have 
higher  resolution  than  systems  with  the  other  laser  sources.  The 
broadband  laser-direct  detection  receiver  system  is  found  to 
have  the  highest  resolution. 

vi 


A system  consisting  of  a quasimonochromatic  laser-heterodyne 
receiver  is  used  for  ranging.  The  range  from  the  system  to  the 
rough  surface  is  found  by  measuring  the  phase  delay  in  the 
reflected  laser  field.  It  is  found  that  the  ranging  performance 
of  a heterodyne  scanning  system  is  not  fundamentally  degraded 
by  speckle  effects. 


AN  ANALYSIS  OF  THE  EFFECTS  OF 


SPECKLE  ON  LASER  SCANNING  SYSTEMS 


I.  Introduction 


A laser  line-scan  imaging  system  is  a system  which  produces 
the  image  of  some  desired  object.  Fundamentally,  the  system 
is  composed  of  a transmitter  and  a receiver.  The  transmitter 
is  a laser  which  illuminates  the  desired  object  surface.  The 
laser  field  is  reflected  or  backscattered  from  the  surface  and 
partially  collected  by  the  receiver.  The  backscattered  field, 
which  is  a function  of  the  reflectivity  of  the  object  surface, 
is  processed  by  the  receiver  to  form  an  image  of  the  object 
surface.  This  image  is  the  desired  output  of  the  receiver. 

The  receiver  contains  optics,  detectors,  filters,  and  other 
components  which  are  needed  to  produce  the  desired  image.  A 
drawback  to  this  system  occurs  because  most  surfaces  are 
rough  when  compared  to  the  wavelength  of  a laser  field.  The 
rough  surface  introduces  noise  into  the  image.  For  coherent 
illumination  this  noise  is  commonly  called  speckle,  and  it 
distorts  the  image  obtained  by  the  receiver.  Speckle  is  discussed 
in  detail  in  Chapter  II. 

The  purpose  of  this  thesis  is  to  determine  the  effects  of 
speckle  on  the  performance  of  a laser  line-scan  imaging  system. 
Three  laser  sources:  monochromatic,  quas imunochromat ic , and 
broadband;  and  two  receivers:  heterodyne  and  direct  detection 
are  investigated.  In  this  thesis,  the  effect  of  propagation 
on  all  fields  is  determined  from  the  Huygcns-Frcsnel  integral. 
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The  rough  surface  is  represented  by  a well  known  statistical 
model  and  the  receiver  output  is  a current  which  is  based 
on  a known  receiver  model.  The  output  current  is  described 
by  a second  moment  mode_  which  is  a function  of  the  laser  beam 
spot  size,  the  system  scanning  velocity,  the  surface  area 
of  the  receiver  input,  and  the  type  of  laser  used. 

The  thesis  is  organized  in  the  following  manner.  Chapter 
II  developes  the  background  material  necessary  for  examining 
the  line-scan  imaging  systems.  In  Chapter  III,  the  second 
moments  models  are  determined  for  two  systems:  monochromatic 
laser-direct  detection  receiver  and  monochromatic  laser- 
heterodyne  receiver.  These  systems  are  compared  on  the  basis 
of  their  signal  to  noise  ratio  and  resolution  ability.  This 
chapter  establishes  the  method  for  determining  the  second 
moment  models  of  the  other  systems.  In  Chapter  IV,  linear 
frequency  modulation  is  used  to  produce  a quasimonochromatic 
and  a broadband  laser  source.  The  results  of  Chapter  III  are 
extended  to  produce  four  more  second  moments  models,  which 
arise  from  the  two  modulated  laser  sources.  The  signal  to 
noise  ratio  and  the  resolution  ability  of  each  model  is 
discussed.  In  Chapter  V,  a quasimonochromatic  laser  source  is 
amplitude  modulated.  The  modulated  laser  field  is  propagated 
to  the  rough  surface  and  back,  and  its  change  in  phase  is 
measured  by  a heterodyne  receiver.  The  phase  change  is  used 
to  determine  the  distance  from  the  surface  to  the  system. 

The  effects  of  quantum  noise  are  also  included  in  this  system. 


1 1 . General  Configuration 


In  this  chapter,  several  background  areas  are  discussed. 
These  areas  serve  as  a basis  for  examining  the  laser  line- 
scan  imaging  system.  The  background  areas  are:  a discussion 
of  speckle,  the  complex  baseband  representation  for  optical 
fields,  the  optical  receiver  models,  the  second  moment  model, 
the  Huygens  - Fresnel  integral,  and  propagation  and  reflection 
of  a laser  field. 

Speckle 

When  a coherent  optical  field  is  incident  on  a surface 
which  is  rough  in  comparison  to  the  optical  wavelength,  a 
random  interference  pattern  is  produced  in  the  reflected  field. 
This  random  interference  pattern  is  called  speckle.  One 
example  of  speckle  is  the  sparkling  pattern  produced  when  the 
light  of  a visible  laser  is  reflected  from  a wall.  The  effect 
of  speckle  is  to  degrade  the  image  of  the  object  surface  by 
superimposing  a "noiselike  structure  which  masks  the  spatial 
information  present  in  the  image"  (Ref  1:1257).  This  results 
in  a pattern  of  constructive  and  destructive  interference 
being  produced  in  the  light  reflected  from  the  rough  surface. 

The  random  interference  pattern  is  composed  of  bright  spots 
due  to  constructive  interference,  dark  spots  due  to  destructive 
interference,  and  areas  with  intensities  between  these  extremes. 
Thus  the  pattern  has  a granular  appearance. 

Rough  surfaces  are  not  the  only  cause  of  speckle.  The 
term  speckle  has  been  generalized  to  include  most  spatial 


interference  effects  that  result  when  any  type  of  wave  or 
field  is  scattered  from  diffuse  objects  (Ref  2:16).  Speckle 
also  occurs  in  atmospheric  transmissions  and  holography.  It 
was  once  thought  to  be  a nuisance,  but  now  it  has  many 
applications.  For  example,  it  is  used  in  the  measurement 
of  motion  of  a rough  body  (Ref  1:1271),  photographic  optical 
processing  (Ref  1:1275),  and  the  remote  sensing  of  crosswind 
in  the  atmosphere  (Ref  3:1). 

The  speckle  phenomena  has  been  known  and  investigated 
since  the  time  of  Newton  (Ref  4).  Its  early  history  is 
briefly  discussed  in  the  introduction  of  a book  edited  by 
Dainty  (Ref  5:1-7).  The  invention  of  the  first  CW  laser 
in  1960  led  to  the  first  observations  of  laser  speckle.  They 
were  reported  by  Rigdcn  and  Gordon  (Ref  6:2367-2368)  and  by 
Oliver  (Ref  7:220).  Since  the  early  1960's,  much  work  has 
been  done  in  the  speckle  related  field.  The  entire  issue  of 
the  November  1976  Journal  of  the  Optical  Society  of  America 
(Ref  1)  was  devoted  to  the  discussion  of  numerous  speckle 
phenomena  and  effects. 

Complex  Baseband  Representation 

Consider  the  real  scalar  optical  field 

u(r,t)  = A(r,t)  cos[2irfot  - <t>(r,t)]  (1) 

It  is  a function  of  time,  t,  and  space,  r,  where  r is  a vector 
with  spatial  coordinates  (x,y,z).  The  complex  baseband 
representation  of  u(r,t)  is 


(2) 


u(r,t)  = Re  { A ( r , t ) exp  [ j <J>  (r  , t)  ] exp  [- j 2tt  f Qt  ] } 

= Re{U(r,t)  exp  [- j 2TrfQt]  } 
where  Re{*}  is  the  real  operator  and  U(r,t)  is  the  complex 
envelope  of  u(r,t). 

Throughout  this  thesis,  the  field,  u(r,t) , is  represented 
by  its  complex  envelope.  The  exponential  time  function  is 
suppressed.  In  addition,  the  field  is  often  described  at  some 
plane  z,  thus  it  varies  in  the  x and  y directions  only. 

Therefore,  the  real  scalar  optical  field  at  a given  plane  z 
is  denoted  U(x,y,t).  The  real  field  can  always  be  obtained 
from  Eq.  (2) . 

Optical  Receivers 

An  optical  receiver  is  a device  which  measures  optical  fields. 
The  receiver  input  is  a field  and  its  output  is  a signal 
proportional  to  the  input  field.  The  receiver  input  for  a line- 
scan  imaging  system  is  the  laser  field  backscattered  from  the 
object  surface.  The  field  entering  the  receiver  is  in  general 
limited  by  the  receiver's  optics,  field-of-view,  or  detector 
area.  It  is  assumed  in  this  thesis  that  the  receiver  is  designed 
so  that  the  field  is  limited  only  by  the  surface  area  of  the 
detector.  This  insures  that  the  entire  image  of  the  object 
surface  and  its  spatial  frequency  content  enters  the  receiver. 

Two  receivers,  the  heterodyne  and  the  direct  detection,  are 
investigated  in  this  paper.  With  the  exception  of  Chapter  V, 
both  receivers  are  assumed  to  be  ideal.  That  is,  all  quantum 
effects  and  noise  terms  are  neglected.  This  allows  the  effects 
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of  speckle  on  the  receiver  output  to  be  isolated.  This  does  not 
make  the  problem  any  less  general.  The  quantum  effects  and 
noise  terms  can  be  added  to  the  final  results  if  desired. 


The  Heterodyne  Receiver.  The  heterodyne  system  adds  the 
input  field  to  a local  oscillator  field  centered  at  the  optical 
frequency  fQ-fjp.  The  ideal  receiver  output  current  is 
CRef  8:481-487) 

iff)  = | j |Us(x,y,t)  + ULQ(x,y)  exp  [j  2irf  ipt]|  2dxdy 


■ f \ Musc*.y.t) 


+ luL0(x,y)l 


+ 2Re(Us(x,y ,t)UL0(x,y)  exp  [- j 2irf Ift] }]dxdy  (3) 


where 

q is  the  charge  of  an  electron 
n is  the  detector  quantum  efficiency 
h is  Planck's  constant 
fQ  is  the  optical  frequency 

is  the  area  of  the  detector  surface 
Us(x,y)  is  the  signal  field  incident  on  the  detector 
surface 

The  first  and  second  terms  in  Oq.  (3)  are  centered  at  zero 
frequency.  They  can  be  electrically  filtered  out  leaving 

ih(t)  = IrP  f f Re{Us(x,y)U*0(x.y)  exP  [*  j 2Tf  Ipt]  Jdxdy  (4) 
o \ 1 


From  Eq.  (4)  it  is  seen  that  the  output  current  is  related 
to  the  actual  input  field  Us(x,y). 

The  Direct  Detection  Receiver.  The  direct  detection 
system  is  a square  law  device.  The  ideal  receiver  output 
current  is  (Ref  9:91) 

id(t)  = t3t"  { | 1 2jxdy  (5) 

0 Ad 

In  Eq.  (5)  it  is  seen  that  the  output  current  is  proportional 
to  the  intensity  of  the  input  field,  |Us(x,y)|  . Recall 
that  the  output  of  the  heterodyne  receiver  is  proportional  to 
the  input  field  itself.  This  is  the  major  difference  in  the 
mechanics  of  the  two  receivers. 

Second  Moment  Model 

The  input  field  to  the  optical  receiver  may  be  unknown. 

This  is  true  for  the  case  where  speckle  occurs  since  the 
receiver  input  field  contains  a random  interference  pattern. 

If  the  input  field  is  unknown  the  receiver  output  current  is 
also  unknown. 

In  statistics,  one  representation  for  an  unknown  is  the 
second  moment  model.  In  this  model  the  unknown  is  thought  of 
as  a signal  plus  an  additive  noise.  The  "signal”  is  considered 
to  be  the  mean  of  the  unknown.  For  an  unknown  current,  the  mean 
is  denoted  E [ i ( t ) ] where  E[*]  is  the  expected  value  operator. 

The  "noise"  is  considered  to  have  a mean  of  zero  and  the  same 
covariance  as  the  unknown.  The  covariance  of  an  unknown 
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current  is  defined  as  (Ref  10:321) 

- E[(i(t)  - E [ i (t)  ] ) (i  (t  ')  - E [ i (t  ")  ] ) ] 

- E [i (t) i (t ") ] - E [i (t) ]E[i  (t ')  ] 

^ R.(t,t')  - E (i  (t)  ] E [i  (t  ■*)  ] (6) 

where  R^(t,t')  is  the  correlation  function  (Ref  11:337).  The 
period  for  which  the  correlation  function  is  non-zero  is  called 
the  correlation  interval.  The  second  moment  model  for  an 
unknown  current  is  illustrated  in  Fig.  1. 

The  Huygens  - Fresnel  Integral 

Consider  the  geometry  of  Fig.  2.  The  complex  envelope  of 
the  monochromatic  field  at  a point  (x,y)  in  the  observation 
plane  due  to  complex  envelope  of  the  field  in  the  object  plane 
is  expressed  by  the  lluygens-Fresnel  integral  (Ref  12:58) 
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Fig.  2 Diffraction  Geometry  (Ref  12:57) 


U2(x,y)  *||  h (x  , y ,a , 3)  ^(a.B)  dad0 
where 

1 exp[jkr01]  

h(x,y,a,B)  » cos(n,rn.) 

j*  r0i 

and 

k is  the  wave  number,  ^ 

Tqj  is  shown  in  Fig.  2 

cos(n,rgj)  is  the  cosine  of  the  angle  between  r^ 
and  the  normal  to  the  a-B  plane. 

Eq. (7)  is  valid  for  monochromatic  fields  only.  The  finite 
extent  of  the  field  in  the  a, 8 plane  is  included  in  the 


mathematical  description  of  U(a,B).  This  restriction  allows 
the  infinite  limits  to  be  placed  on  the  integral  in  Eq.  (7). 

Unless  otherwise  noted,  all  integrals  in  this  paper  have 
infinite  limits.  For  simplicity  these  limits  are  suppressed. 

Eq.  (7)  is  difficult  to  use.  It  can  be  simplified  through 
several  approximations.  If  the  angle  between  n and  r^ 
is  less  than  18°,  then  cos(n,r01)  = 1.  Also,  if  the  deviation 
of  the  field  from  the  z axis  is  much  less  than  the  distance 
z,  the  quantity  r^^  in  the  denominator  of  Eq.  (8)  is  approximately 
z.  Eq.  (8)  now  becomes 

h(x,y,o,B)  = exp[jkr01]  (9) 

The  r^  in  the  exponent  can  not  be  approximated  by  z 
because  it  is  multiplied  by  a large  number,  k.  Thus  even 
small  changes  in  can  result  in  phase  changes  much  greater 
than  2tt. 

One  commonly  used  exponential  approximation  is  the  Fresnel 
approximation.  From  Fig.  2 

roi  = + ( + ( YrL)2}h  (10) 

In  the  Fresnel  approximation,  the  binomial  expansion, 

[l+b)h  = 1 + \ b - | b2  + •••  | b | < | (11) 

is  used.  Only  the  first  two  terms  in  the  expansion  are  kept. 

Thus  the  exponential  approximation  for  r^  is 

r01  ■ 1 t1  * I<  ^)2  * r C (12) 
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Goodman  (Ref  12:59)  states  that  this  approximation  is  good  when 


z > > -jy  [(x-a)  + (y-B)  ], 


Now  Eq.  (9)  can  be  written  as 


h(x,y,a,S)  = ^-p2k2]  exp[i|[(x-a)2  + (y-g)2]]  C 

Eq.  (7)  then  becomes  the  Fresnel  approximation  written  as 
U (X  y)  = eXp  [ j kz]  f 

U2tx*yJ  JAz  J U^a.g) 


exp[^|  [(x-a)2  + (y-B)2]]  dadg 


U2(x,y)  - ~ exp  [j-^-(x2+y2)  ] j j U^a.g) 

exp[j-y£-  (a2  + g2)  ] exp  [- j-£-(xa  + yg)]  dadg  (16) 


In  many  cases,  field  propagation  is  adequately  described  by 
the  one  dimensional  form  of  Eq.  (16).  It  is  given  by 
(Ref  13:316) 

exp[j  (kz-7^  )]  . 2 

U,(x)  = T - exp[j  ] U (a) 

(Az) 2 J J 1 

exp[j^-]  exp[-j^|^]  dx  (17) 

In  the  latter  case,  the  field  at  the  z plane  is  variable  in  only 
one  lateral  direction. 

As  shown  in  Appendix  A,  the  IIuygens-Fresnel  integral  can 
be  extended  to  include  the  propagation  of  broadband  fields. 

The  one  dimensional  extended  integral  is 
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it  2 

U2(x,t)  = exp  [ j (kz  - T )]  exp [ j k 

(-*,)■* 

T 2 

| | Uj  (a , t - -£-)  exp[j^]  exp[-j2TTj|]da  (18) 

r01 

where  c is  the  speed  of  light  and  the  quantity  /c  is  due 
to  the  propagation  delay  of  the  field.  Eq.  (18)  is  valid 
for  all  time  varying  fields  as  long  as 


o 


(19) 


where  B is  the  bandwidth  of  the  complex  envelope  which  is 
centered  at  the  frequency  fQ. 

A special  case  of  Eq.  (18)  is  valid  for  the  propagation 
of  quasimonochromatic  fields.  The  bandwidth  restriction  on 
these  fields  is  that  (Ref  12:108) 


1 > J>  | r0 1 1 max 


(20) 


This  restriction  implies  that  for  a fixed  a,  U(a,t)  cannot 
change  significantly  in  any  time  much  less  than  seconds. 
Eq.  (18)  is  modified  to  represent  the  propagation  of 
quasimonochromatic  fields  through  the  following  substitution 
(Ref  12  : 5 S ) 


Ujtoi.t 


= Uj(a,t) 


(21) 


Propagation  and  Reflection  of  a Laser  Field 

Siegman  has  shown  (Ref  14:306)  that  when  the  output 
of  a monochromatic  laser  is  propagated  by  use  of  the  Muygcns- 
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Fresnel  Integral,  the  resulting  laser  field  is  a Gaussian  spherical 
wave  of  the  form 

2 2 

U(a,B)  = Aexp[jkz]exp[jj^(a2  + 62)]exp[-  - — 1-1 — ] (22) 

LZ  W*(z) 

where  A is  the  amplitude  of  the  field  and  W(z)  is  the  radius 
of  the  beam  spot  size  at  a distance  z from  the  laser.  The 
radius  is  measured  at  the  point  where  the  field  amplitude  is 
e ^ of  it  maximum  value.  The  one  dimensional  form  of  Eq.  (22) 
is 

2 

Ufa)  = explj  (kz-JOjexp  J— £ ]exp[j  a2]  (23) 

(Xz)J  * Vv  (z)  lz 

When  the  above  field  is  normally  incident  on  a surface  that 

is  smooth  with  respect  to  X,  it  is  reflected  by  conjugating 

the  exponential  terms  and  multiplying  the  resulting  field  by 

the  surface's  reflection  coefficient.  When  exp[jkz]  is  conjugated, 

it  implies  that  the  field  is  progating  in  a direction  which  is 

opposite  to  that  of  the  incident  field.  Conjugating 
k 2 

exptj-^y-  x ] inverts  the  spherical  wavefront  of  the  field. 

Some  of  the  incident  field  may  be  absorbed  by  or  transmitted 
through  the  surface.  To  represent  this  loss,  a surface  reflection 
coefficient  is  used.  This  coefficient  has  a value  between 
zero  (no  field  reflected)  and  one  (entire  field  reflected) . 

The  reflected  field  is  written 

Ur(a)  = a(a)U.*  (a)  (24) 

where  a(ct)  is  the  surface  reflection  coefficient  and  the 
subscripts  r and  i indicate  the  reflected  and  incident  fields, 
respectively. 
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The  development  of  the  background  material  is  now  complete. 
The  ideas  and  theory  discussed  here  are  used  to  examine  the 
line  scan  imaging  system. 


III.  The  System  Models  for  a Monochromatic 
Laser  Source 

The  method  for  determining  the  system  models  is  developed 
in  this  chapter.  First  a representation  for  the  mechanics  of 
the  imaging  system  is  discussed  and  the  rough  surface  model  is 
defined.  These  ideas  are.  combined  with  the  background  to 
determine  the  system  models  for  a monochromatic  laser- 
heterodyne  receiver  system  and  a monochromatic  laser-direct 
detection  receiver  system.  For  simplicity,  the  models  are 
found  for  one  dimension.  Two  dimension  models  are  straightforward 
but  tedious. 

The  Line  Scan  Imaging  System  Model 

In  a line  scan  imaging  system,  the  laser  scans  across  some 
object  surface.  The  mechanics  of  the  system  keep  the  detector's 
surface  area  aligned  normal  to  the  direction  of  the  laser  field. 

As  the  object  surface  is  scanned,  the  angle  between  the  incident 
laser  field  and  the  surface  changes.  The  analysis  of  the  system 
is  simplified  if  the  system  mechanics  and  change  in  incident 
angle  are  ignored.  These  two  conditions  are  satisfied  under 
the  following  two  assumptions:  (1)  the  laser  field  is  normal 
to  the  object  surface  at  all  times,  and  (2)  the  propagation 
of  the  lasher  field  is  normal  to  the  detector's  surface  at  all 
times . 

During  typical  use,  the  line  scan  imaging  system  is  mounted 
in  an  airplane  which  flies  over  the  object  surface.  To  represent 
the  movement. of  the  system  past  the  surface  it  is  assumed  that 
the  coordinate  system  of  the  laser  and  receiver  is  fixed  in  space, 


Detector 

Surface 


a (a+vt)' 


Fig.  3 Line  Scan  Imaging  System  Model  (Ref  16:28) 


while  the  coordinate  system  of  the  surface  and  its  reflection 
coef f icient , ' a (x) , is  traveling  in  time  with  velocity  v.  The 
model  is  a function  of  the  velocity  and  the  laser  beam  width, 
W(z),  at  the  object  surface.  By  varying  these  two  parameters 
it  is  possible  to  determine  the  effects  of  spreading  beam 
width  and  changes  in  surface  velocity.  The  system  model  is 
illustrated  in  Fig.  3. 


The  Rough  Surface  Model 


The  field  reflected  from  a rough  surface  is,  in  general, 
difficult  to  determine.  It  is  dependent  on  variations  in 
surface  height,  the  polarization  of  the  incident  field,  and 
the  reflection  coefficient  of  the  surface.  In  order  to  avoid 


a detailed  study  of  the  surface  scattering  problem,  a simple 


model  is  developed  to  represent  the  rough  surface  (For  a deeper 
mathematical  understanding  of  the  properties  of  electromagnetic 
fields  scattered  from  a rough  surface  consult  Beckman  and 
Spizzichino  (Ref  17)).  'the  model  is  developed  under  five 
conditions . 

First,  depolarization  effects  are  ignored  becasue  the 
polarization  of  the  incident  field  is  not  included  in  this 
thesis.  Second,  multiple  scattering  at  the  surface  is 
insignificant  since  only  the  field  reflected  normal  to  the  object 
surface  is  incident  on  the  detector.  These  two  conditions  are 
typically  ignored  (Ref  17:  Chapters  1,  3,  and  5;  18:1689) 
since  their  effects  are  very  small  and  difficult  to  describe 
mathematically.  Third,  the  surface  is  assumed  to  be  rough 
when  compared  to  the  optical  wavelength  so  that  speckle  effects 
occur.  The  fourth  condition  is  that  the  rough  surface  is 
thought  of  as  a random  variation  in  height  (Ref  17:  Chapter  3). 

This  allows  the  surface  roughness  to  be  described  by  statistics. 
Also,  a field  incident  on  this  surface  undergoes  a random  variation 
in  phase.  The  reflected  field  is  then  equivalent  to  a field 
reflected  from  a smooth  surface  and  multiplied  by  a random  phase 
term.  Finally,  the  surface's  reflection  coefficient  is  assumed 
to  be  a real  function  of  space  with  a value,  between  zero  and 
one  so  that  it  attenuates  the  field.  The  reflection  coefficient, 
a(a),  is  the  quantity  that  is  measured  by  the  system.  These 
conditions  are  the  same  as  the  ones  used  by  Lyons  in  a similar 
heterodyne  imaging  problem  (Ref  16:21).  The  model  is  also 
similar  to  the  ones  used  by  Miller,  et.  al.  (Ref  19:779-785) 
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and  Beckman  and  Spizzichino  (Ref  17:  Chapter  5). 

Recall  that  Eq.  (24)  describes  a monochromatic  field 
reflected  from  a smooth  surface.  The  field  reflected  from  a 
rough  surface  is  determined  by  applying  the  five  rough  surface 
model  conditions  to  Eq.  (24).  The  reflected  field  is 

UrCa)  « a (a)  exp  [j  0(a)]  u/(a)  (25) 

where  0(a)  is  the  random  phase  and  is  a function  of  space. 

The  relation  between  the  random  variations  in  surface  height, 
H(x) , and  the  random  phase  was  previously  shown  to  be 
(Ref  20:36;  21:157) 

0ta)  = H(a)  (26) 

The  variation  in  surface  height  is  modeled  as  a zero  mean, 
stationary,  Gaussian  random  process.  The  zero  mean  implies  the 
point  that  the  heights  are  measured  from  is  chosen  in  a manner 
such  that  E[H(a)]  = 0.  The  rough  surface  is  often  described  by 
a Gaussian  distribution  (Ref  5:65;  1:1153,  1195,  1205,  1212, 
1224;  17:80).  Using  these  assumptions,  it  is  possible  to 
determine  the  mean  and  variance  of  the  phase.  The  mean  is 

E[0(o)]  = il  E [H  (a)  ] = 0 (27) 


and  the  variance  is 


2 

ae  = 
where 

roughness 


C^)2  E [H2  (a)  ] = Cx^aH2  (28) 

is  the  variance  of  the  surface  height.  The  RMS  surface 
a^,  can  be  compared  to  the  optical  wavelength,  A , 


by  use  of  Eq.  (28).  If  the  surface  roughness  is  small  compared 
to  the  optical  wavelength,  the  RMS  phase  variation,  oQ , is 
also  small.  This  means  there  is  little  problem  with  interference. 
However,  as  the  surface  roughness  approaches  the  optical  wave- 
length, the  RMS  phase  variation  approaches  2tt.  Thus  interference 
is  produced  in  the  reflected  field.  This  interference  is 
called  speckle. 

The  Direct  Detection  Current  Model 

The  ideal  detector  output  current  for  a direct  detection 
system  can  now  be  determined.  Recall  that  this  current  is  due 
to  the  field  input  to  the  detector.  The  one  dimensional 
monochromatic  laser  field  incident  on  the  rough  surface  is 
given  by  Eq.  (23) 

a 2 . 2 

U.  (a)  = — A— p- - exp  [ j (kz- j)  ] exp  [- — § ] exp[jJ~  ] (29) 

1 (Xz) a 4 tr(z) 

The  field  reflected  from  the  rough  surface  is 

Ur(cO  = a (a  + vt)  e^ 0 (a+vt)  (a)  (30) 

The  argument  "ot+vt"  results  becasue  the  surface  and  thus  its 
reflection  coefficient  and  height  variation  is  moving  with 
respect  to  the  laser.  Substituting  Eqs . (29)  and  (30)  into 
Eq.  (17)  gives  the  field  at  the  input  to  the  detector.  It  is 
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UdCx)  = Tz~  exPt*j^Ckz-7T/^)] exp [ - j k 

f 2 

a (a+vt)  exp  [ j 0 (a+vt)  ] exp[-  — ] 
J wz(x) 


exP  [*  j^~ ] exp  [ j 2tt^]  da 
The  magnitude  squared  of  this  field  is 


(31) 


|U  , Cx)  | 


<TT>2 


2 + ft2 

a (a+vt) a (B+vt) exp  [ - ^ ] 

" (z) 


2 2 

exp[-jk  — ] exp  [ j 2n-^~ - (a-B)] 


exp [j 0 (a+vt) - i 0 ( B + vt) ] dad3 


(32) 


Eq.  (32)  represents  the  intensity  of  the  field  at  the  detector 
input.  The  ideal  receiver  output  current  for  the  direct  detection 
system  is  given  in  Eq.  (5).  Its  one  dimensional  form  is 


idCt)  = f V*)|Us(x)|2dx 


(33) 


where  PD(x)  is  the  limiting  aperture  function  which  describes 
the  surface  of  the  detector.  Combining  Eqs . (32)  and  (33) 
yields  the  detector  output  current  which  is 


V1’  ■ <-&-> 2 h¥ 


ia 

• 

r 

K 

i 

PD(x)a(a+vt)a(B+vt) 


2 a1  2 2 

exp[-  *2*  ^ ]exp[-jk°*  ]exp[j^^(a-g)] 


W (x) 

exp[j (0(a+vt) -e(0+vt))]  dxdadB 


(34) 
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To  simplify  the  notation  let 


and 

2 

h(a)  = exp  [ - -J5 — ] 

• W^U) 


(35 -a) 


(35 -b) 


where  h(a)  is  called  the  system  function.  It  determines  the 
effect  of  the  beam  spot  size  on  the  system.  Solving  the  x 
integration  in  Eq.  (31)  yields  the  spatial  Fourier  transform 
of  P (x)  evaluated  at  f^  = . The  transform  of  P^(x) 

will  be  denoted 


Fx[pD(x)] 


q-g 

Xz 


^DF  O 


(36) 


Eq.  (34)  now  reduces  to 

id  CO  = A1  ||  PDF(a-g)  a(a+vt)a(8+vt) 

2 2 

h(a)h(8)  exp[-jk^-ll — ] 

exp[j (6 (a+vt) -0 (B+vt) ) ] dadB  (37) 

The  Mean  Detector  Current.  Now  that  the  detector  current 
has  been  determined,  it  is  possible  to  find  the  mean  or  "signal" 
component  of  the  current.  The  only  random  term  in  Eq.  (37) 
is  the  phase  deviation,  thus  the  mean  of  the  current  is 
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E[id(t)]  = Ax  j PDF(a-S)a(a+vt)a(B+vt) 


2 p2 

h(a)  h ( 6 ) exp  [ - j ka-  -~-b — ] 


E [exp (j  (9 (a+vt)  -0(3+vt))]  dadg 


(38) 


The  phase  deviation,  0,  is  a zero  mean  Gaussian  random  process 
Therefore  the  term,  E [exp [j (0 (a+vt) •- 0 (B+vt) ] ] can  readily  be 
determined  by  use  of  the  characteristic  function.  A form  of 
the  characteristic  function  (Ref  10:419)  is 


(V)  = E[ejVx] 


(39) 


If  x is  a Gaussian  random  variable,  then  Eq.  (39)  becomes 

(40) 


v2o  2 

<(>x(V)  = exp  [ j VE  [x]  - — ] 


In  Eq.  (39)  let  V = 1 and  x = 0(a+vt)  -0(B+vt).  Then 


2_  2- 


E [exp  [j(0(a+vt)-0(  B+vt)  ) ] = exp[-*s(l)  ] 


where 


a12  = E [ (0 (a+vt)  - 0 ( 8+vt ) ) 2 ] 

= E [02 (a+vt)  + 02(0+vt)  -20 (a+vt) 0 (B+vt) 


2 CSr 


2Rg (a- 6) 


£ 


2 o0  (1- p (a- 8) 


(41) 


(42) 


where  p(a-B)  is  the  normalized  correlation  function  of  the 
phase  0.  Eq.  (41)  now  becomes 
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1 


1 


E[exp[j (6(a+vt)  -6(B+vt))]  = exp[-aQ  (l-p(a-B))]  (43) 

This  function  has  a value  of  one  when  oc-0=  0,  and  decreases 
2 

to  exp[-Og  ] as  ot-B  approaches  infinity. 

Beckman  and  Spizzichino  (Ref  17:81)  suggest  that  the 
phase  correlation  function  can  be  represented  by  a Gaussian 
function  of  the  form 


p(a-B)  = exp  [ -(^)2] 

c 


(44) 


where  r is  the  correlation  distance  of  the  rough  surface 
when  p(ct-B)  drops  to  the  value  e Goodman  (Ref  18:1698) 
states  that  rc  is  generally  less  than  0.1mm  for  optically 
rough  surfaces. 

It  has  been  assumed  that  the  deviation  of  the  surface 
roughness  is  greater  that  the  optical  wavelength  or  o^  > X. 
Eq.  (28)  can  be  rewritten  as 


4lT  ^ 4 TT  , . 

°e  " X > T x ' 4,1 


(45) 


When  crQ  > 5,  the  correlation  distance  of  the  reflected  field 

is  less  than  the  correlation  distance  of  the  rough  surface 

by  a factor  of  1/  (Ref  22:984).  Therefore,  the  field 

9 

correlation  distance,  l,  is 


A < T°/, 


1mm  - n i 

= . 01mm 


'0 


4? 


(46) 


This  field  correlation  distance  was  previously  derived  by 
Lyons  (Ref  16:33).  Its  significance  is  that  it  is  the  maximum 


l 

y 
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distance  that  any  two  points  in  the  reflected  field  (at  the 
surface)  can  be  separated  by  and  still  be  correlated.  Eq.  (43) 
is  now  called  the  field  correlation  function  and  is  denoted 
(a- 3) . That  is 

P^fot-B)  = exp [-Og2(l-p(a-$))]  (47) 

where  J i indicates  the  field  correlation  distance.  Fig.  4 

Ot  ~ R 

is  a plot  of  P.(a-B)  as  a function  of  — - for  several  values 

rc 

of  a^j. 

Another  simpliciation  to  Eq.  (38)  is  the  approximation 

exp [- j|(a2- B2) ] = 1 (48) 

which  is  valid  for  the  condition  of 


— | ot2  - 62  | < < 1 

z1  'max 


(49) 


2 2 

|a  -B  I occurs  when  a takes  on  its  maximum  value  and  6 
1 'max 

is  at  its  minimum,  or  vice  versa.  The  maximum  value  of  a 
and  B is  the  radius  of  the  beam  spot  size,  W(z) , on  the  rough 
surface.  The  maximum  value  that  a and  B can  be  separated 
and  still  be  correlated  (i.e.  nonzero)  is  l,  the  field 
correlation  distance.  Thus 


I a2 -6  2 | max  = UW(z))2  ‘ (W(z) -£)  2 | 
= 2W(z)JL 


(50) 


Eq.  (49)  becomes 

£ 2W(z)  8.  < < 1 

z 


(51) 
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7 


aH  = X 


Fig.  4 Plot  of  P£(a-B)  Ref  (16:34) 


which  can  be  rewritten  as 


WCH  < < * 

2 4tt£. 


„ = X 

aH  ' & 


where  W-z--  is  the  beam  divergence  of  the  laser.  The  laser 
can  be  designed  to  meet  the  condition  of  Eq.  (52).  For  example 
the  beam  divergence  of  a 1.06ym  laser  must  be  less  than  8 
milliradians . Beam  divergences  are  typically  less  than  1 
milliradian,  thus  the  approximation  of  Eq.  (48)  holds.  Now 


the  mean  of  the  current  becomes 

E[id(t)]  - Ax  | | PDF(a-8)a(a+vt)aCe+vt) 

h(a)h(8)Pa (a-8)dadB  (53) 

Due  to  the  approximation  in  Eq.  (48),  the  mean  of  the  current 
is  no  longer  dependent  on  the  quadratic  terms  in  the  exponential. 
Thus,  the  resulting  mean  is  identical  to  that  obtained  in  the  far 
field  case. 

In  Eq.  (53),  only  the  field  correlation  function,  P^(a-8) 
and  the  spatial  Fourier  transform  of  the  detector  surface, 
Ppp(a-8),  are  dependent  on  the  argument  "a-B".  Recall  that 
P^(a-B)  was  plotted  in  Fig.  4.  For  the  case  of  = A,  P^fot-B) 
is  plotted  once  more  in  Fig.  5a.  Note  that  the  width  of  P^a-B) 
is  only  a few  micrometers.  Generally,  the  width  Ppp(a-B)  is 
much  wider.  For  example,  assume  that  the  detector  surface  is  of 
the  form  in  Fig.  5b.  It  has  a width  D,  and  a uniform 
transmittance  of  one.  Then 


PDF<°-«  = 


■f  = a~  B 
Xx  - Az 


sin(27r«^  u) 

ajl 

Az 


(54) 


PDF(a-B)  is  plotted  in  Fig.  5c.  The  width  of  Ppp(a-B) (ignoring 

X z 

the  side  lobes)  is  /q/Z'  For  tyP*031  values  of  A = 1.06ym, 
z = 103  meters,  and  D = 10  centimeters,  the  width  of  PDp(a-8)  is 
approximately  20  millimeters.  Because  PDp(a-  8 ) is  much  wider  than 
P^fa-B),  it  is  constant  over  the  region  where  P^(a-B)  is  non-zero. 
Thus  their  effective  product  is 


pDfC«-HV“-«  = PDF(0)Pt(a'6) 


(55) 
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Eq,  (53)  now  becomes 

E [id Ct)  ] = A^DpCO)  | [ a (a+vt) a (B+vt) 

h(a)h(B)  P^(a-B)dadB 


(56) 


Using  the  change  of  variables,  x = a + vt  and  y = B+vt, 
the  mean  is 


E[id(t)] 


aipdf^°-) 


a (x)  a (y)  P 


(x-y) 


h(x-vt)h(y-vt)  dxdy 


(57) 


Eq.  (57)  is  the  "signal"  component  in  the  second  moment 
model  for  a direct  detection  receiver.  The  model  will  be 
presented  after  the  covariance  of  the  receiver  current, 

Eq.  (37),  is  found.  Eq.  (57)  is  valid  for  the  following 
conditions:  (1)  the  field  incident  on  the  detector  is  due  to 
the  reflection  of  a monochromatic  laser  field  from  a rough  surface, 

(2)  the  RMS  surface  roughness  is  greater  than  the  laser  wavelength, 

(3)  the  laser  is  design  so  that  its  beam  divergence  is  less 
than  j—  , and  (4)  PDp(ct-B)  is  much  wider  than  P^(a-B).  Now 
the  covariance  will  be  found. 

The  Covariance  of  the  Current.  In  the  second  moment  model, 
the  covariance  describes  the  "noise"  or  the  fluctuation  of  the 
current  from  its  mean.  It  is  found  by  first  determining  the 
correlation  of  the  current.  From  Eqs.  (6)  and  (37) 
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R,  Ct.t-) 

*d 


E[id(t) id(t ') ] 

•ill 


PDF(a-e) PDFCa'-B')a(a+vt) 


a (0+vt) a (a '+vt  *)a(p'  + vt')h(a)h(B)h(a')h(B>) 

E [exp [ j (9 (a+vt)  + 9(a'+vt')  - 0 (B+vt) - 0 (B '+vt ') ) ] 

dadBda'dB' 


] 

(59) 


The  approximation  of  Eq.  (48)  was  used  in  Eq.  (59). 

The  expected  value  of  the  phase  deviation  term  in  Eq.  (59) 
is  determined  through  the  application  of  the  moment  factoring 
theorem.  Reed  (Ref  23:  194-195)  has  shown  that  if  Z^,  Z^,  Z^, 
and  Z4  are  zero  mean,  complex  jointly  Gaussian  random  processes, 
then 


e[z1z2z3*z4*j  = e[z1z3* 

+ E[ZFZ4* 


]E[z2z4  ] 

] E[Z2Z3*] 


(60) 


One  problem  results  when  the  moment  theorem  is  applied. 

The  zn  term  in  the  moment  theorem  corresponds  to  a exp[ja] 
term  in  Eq.  (59).  It  was  argued  in  the  section  on  the  rough 
surface  model,  that  9 is  described  by  a Gaussian  distribution. 
Therefore,  the  distribution  of  exp[j9]  is  not  Gaussian.  However, 
the  reflected  field  is  just  a complex  sum  of  the  fields  reflected 
from  each  point  on  the  rough  surface.  It  can  be  shown  by  the 
Central  Limit  Theorem  (Ref  17:191)  that  the  reflected  field 
incident  on  the  receiver  has  a Gaussian  distribution.  Thus, 
it  is  assumed  that  exp[je]  takes  on  a Guassian  distribution 
as  it  propagates  through  the  atmosphere. 
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Therefore,  the  expected  value  term  in  Eq.  (59)  becomes 


<J>0(1)  = E [exp [j (0 (a+vt)  +0  (a  ' + vt  ')  - 0 ( 3+vt)  - 0 (6  '+vt  ') ) ] ] 

*=  E [exp  [j  0 (a+vt)  - j 0 (3+vt)  ] ] E [exp  [j0(a'+vt")-j0  (B'+vt ')  ] ] 

+ E [exp [ j 0 (a+vt) - j 0 (8 >+vt ') ] ] E [exp [ j 0 (a '+vt ') - j 0 (B+vt) ] ] (61) 


The  characteristic  function  can  be  used  to  show  that 
<*>0  (1)  = exp[-a02(l-p(a-B))]exp[-a02(l-p(a'-B'))] 


+ exp [-aQ2 (1- p(a- 6'+  vAt))]  exp [ -afl2 (1 - p (a'- g-vAt ) ) ] (62) 


0 


Recall  that  the  field  correlation  function  was  defined  in 
Eq.  (47)  as 


P£(ot)  = exp  [-Og  (1-P(a))] 


(63) 


thus 


>0(1)  = P£(a-B)P£(a'-0')  + P£(a-B'+vAt)  P£(a'-S -vAt)  (64) 


Recall  further  from  Eq.  (5)  that 

R.  (t,t')  = E [i  , (t)  ] E [i,  (t  ')  ] + C,  (t,t') 


(65) 


The  first  term,  P£ (a- 6) P£ (a 3 ')  , in  the  summation  of  Eq.  (64) 
corresponds  to  E [i^(t) ]E  [i^(t ') ] . Thus  the  covariance  of  the 
current  is 


30 


v 


cid(t,t')  -A l l ' 


r 

a(a+vt)a(a'+vt')a(B+vt) a (3  >+vt ') 


h(a)h(B)h(a')h(B')  PnP(a-B)  Pn_(a'-B') 


P& (a- B '+vAt)  P^ (a'- B-vAt)  dadB  (66) 

Eq.  (66)  is  difficult  to.  interpret  because  the  Fourier  transform 
of  the  aperture  is  dependent  on  "a -B " and  "a'-B'M  but  the  field 
correlation  function  is  dependent  on  "a*B""  and  "a'-6".  Thus  in 
Eq.  (66)  the  following  simplification  is  made: 


^DF^a"^  — PDF^ 


PDF^a  ^ — PDF  ^ ^ (68 

This  simplification  produces  an  upper  bound  on  the  covarince. 
It  is  writteh  .as 


C.  (t,t')  <A? 
d 1 


a(a+vt)  a (B+vt)  a (a  '+vt  ')  a(B  '+vt  ') 


h(a)h(B)h(a')h(6>)  PDF  (0) 


P^(a-B '+vAt)  P^(a'-B-vAt)  dadBda'dB 


which  is  rewritten 


ci  C ^ » t ' ) _<  [A1Pdf(0)  a(a+vt)a(B'+vt')h(a) 


h(B')  P0(a-B'+vAt)  dadB'] 


The  correlation  distance  of  the  current  is  difficult  to 
determine  since  the  covariance,  Eq.  (70),  is  nonstationary. 
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An  idea  of  the  correlation  distance  can  be  obtained  by  writting 
Eq.  (70)  as 

Cid(t,t')  < [AjP^CO)  | a (a+vt) h (a) 

[a (a+vt ") h (a)  * P^ (a+vAt) ] da] ^ (71) 

where  * denotes  the  convoltuion  process.  The  convolution  is 
graphically  represented  in  Fig.  6.  The  maximum  width  of 
P^(a+vAt)  is  l , the  field  correlation  distance.  The  maximum 
width  of  a(a+vt)h(a)  and  a(a+vt')h(a)  is  determined  by  the  beam 
spot  size.  This  distance  is  denoted  as  L.  L is  much  greater 
then  l . When  vAt  is  greater  than  ^ f 2 » convolution  of 

a(a+vt')h(a)  with  P^(x+vAt)  -is  equal  to  zero.  Therefore,  the 
maximum  correlation  distance  of  the  current  is 


■’Ue 


L = 2 vAt  (72) 

and  the  correlation  time  of  the  current  is 


At  = (73) 

From  Eqs . (72)  and  (73)  it  is  seen  that  the  correlation 
distance  increases  with  increasing  correlation  time.  It 
should  be  noted,  that  if  the  reflectance,  a(a+vt'), 
fluctuates  very  much  within  the  distance  L,  then  the  correlation 
distance  and  time  of  Eqs.  (72)  and  (73)  will  become  less. 


The  Direct  Detection  Model.  The  direct  detection  current 
is  modeled  using  the  second  moment  method.  The  model  is 
described  by  the  mean  and  covariance  of  the  current.  The  mean 
was  given  in  Eq.  (57)  as 


E[id(t)]=  Ax  PDF (0) 


a(a)a(P)h(a-vt)h(e-vt) 


PJl(a-8)  dadB 


(57) 


It  was  previously  noted  that  P£(a-8)  is  very  narrow  with 
respect  to  the  other  terms  in  the  integral  in  Eq . (57). 

Due  to  its  narrowness,  P£(a-P)  is  approximated  by  the  Dirac 
Delta  function, 

PA(o-0)  = A26(a-8)  (74) 

where  A 2 is  the  area  of  P£(a-3)  and  is  give  by 


A 


2 


PJl(a-6)d(a-8) 


(75) 


Eq.  (57)  becomes 
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E[id(t)]  = Aj  A2PDp(0)  [ 


a(a)a(B)h(a-vt) 


h(B-Vt) 5 (a  - B ) dads 


Wdf^ 


a2(a)h2(a-vt)da 


(76) 


Recall  from  Eq.  (35)  that  h(a)  is  an  even  function, 
thus , 

E[id(t))  = A1A2Pdf(0)  | a2 (a) h2 (vt - a) da 


A1A2PDF(0) ta2Cx)*h2(x)] 


(77) 


Eq.  (77)  implies  that  the  "signal"  component  of  the  second 

moment  model  is  proportional  to  the  convolution  of  the  surface 

2 

reflection  coefficient  squared,  a (x) , with  the  system  function 
2 

squared,  h (x) . 

To  complete  the  second  moment  model  for  the  direct 
detection  current,  the  "noise"  component  is  found.  Eq.  (70) 
gave  the  upper  bound  on  the  covariance  as 


C £ (t,t')  < [A^pfO)  j [ a(a)a(0)h(a-vt)h(0-vt  ') 
d 


P.  (a-B)dadB] 


(78) 


I 
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This  equation  describes  the  "noise"  term  in  the  second  moment 
model.  The  double  integral  in  the  above  equation  is  equivalent 
to  the  correlation  function  at  the  output  of  a linear  filter 
whose  input  was  a zero  mean  noise.  To  see  this  consider  a 
zero  mean  noise,  n(x),  input  to  a linear  filter,  h(x) . The 
filter  output  is 

nx  (x)  = n(x)*h(x) 

= | n(a)h(x-a)d  (79) 

The  correlation  of  the  filter  output  is 
Rn^(x,x')  = E[n1(x)n1(x")] 

= E [n(a)n(3)  ]h(x-a)h(x'-8)dad6 

. 

= Rn(af3) h (x-a) h (x 3)  dad8  (80) 

Comparing  Eq.  (80)  with  the  double  integral  in  Eq.  (78)  it  is 
seen  the  correlation  function  of  the  noise  is 

Rn ( x ’ x ) = a(x)a(x")Pjl(x-x>)  (81) 

The  square  of  the  integral  in  Eq.  (78)  implies  that  the  covariance 
of  the  current  is  equivalent  the  multiplication  of  two  noise 
sources.  These  sources  are  identically  distributed,  statistically 
independent,  and  have  a correlation  function  given  by  Eq.  (81). 

The  second  moment  model  of  the  direct  detection  current  is  now 
complete.  The  model  is  diagramed  in  Fig.  7. 


Fig.  7 Direct  Detection  Current  Model 


Recall  that  all  noise  terras,  other  than  speckle,  were 
neglected.  Thus  the  noise  in  this  model  represents  speckle 
noise.  The  model  presented  above  can  easily  be  modified  to 
include  the  quantum  noise  effects  that  were  neglected  in  the 
receiver  output  current. 

A parameter  that  is  typically  used  to  measure  performance 
is  the  signal  to  noise  ratio,  denoted  SNR.  The  SNR  is  defined 
to  be  the  ratio  of  signal  power  to  noise  power  in  a system. 

In  statistical  terms,  the  SNR  is  the  mean  squared  divided  by 
the  variance  (Ref  24:264).  For  the  direct  detection  current 
model,  the  SNR  is 
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E2[i(t) 


I 


SNR.  = 
1 


. E2U(t)] 

C(  (82) 

1d 


Substituting  Eq.  (57)  and  (70)  into  Eq.  (82)  gives  the 
result 


SNR.  > 1 (83) 

xd  “ 

The  "greater  than  or  equal  to"  sign  occurs  because  of  the  upper 
bound  on  the  covariance. 

Several  points  can  be  noted  about  the  direct  detection  current 

2 

model.  The  spatial  filtering  effects  of  h (x)  are  seen  by  taking 
the  spatial  Fourier  transform  of  E[i£j(t)]  as  given  in  Eq.  (77). 

Fx[E[id(t)]]  " A1A2PDF(0)A'(fx)[H-(fx)*H'(fx)]  (84) 

where  A'(fx)  - Fx[aZ(x)]  and  H'(fx)  - Fx[h(x)]. 

The  Fourier  transform  has  an  inverse  effect  on  the  width  of  a 

function.  That  is,  as  h(x)  gets  wider  H'(f^)  and  thus  the 

convolution  of  H'(f  ) with  itself  gets  narrower.  The  results  is 

that  as  the  system  function  gets  wider  in  space,  it  blocks  out 

2 

the  high  spatial  frequency  content  of  a (x) . 

The  model  for  E[id(t)]  is  identical  to  the  result  discussed 
by  Goodman  for  an  incoherent  imaging  system  (Ref  13:109). 


Thus  the  direct  detection  current  model  can  be  thought  of  as  an 

incoherent  imaging  system  with  an  additive  noise  term. 

Last  of  all,  it  can  be  seen  from  either  Eq.  (78)  or  Fig.  7 

2 

that  the  noise  is  signal  dependent.  Thus  as  the  signal,  a (x) , 
increases  or  decreases,  the  noise  also  increases  or  decreases  by 
a proportional  amount.  Thus  the  SNR  can  not  be  improved  by 
increasing  the  signal  power.  This  is  in  contrast  to  many  other 
noise  models  which  are  signal  independent. 

This  concludes  the  discussion  of  the  direct  detection  current 
model.  The  current  model  for  the  heterodyne  receiver  is 
determined  in  the  following  section. 

The  Heterodyne  Receiver  Current  Model 

A second  moment  model  for  the  receiver  current  of  a 

monochromatic  laser-heterodyne  receiver  system  has  been 

determined  by  Lyons  (Ref  16:46).  The  procedure  used  to  establish 

the  heterodyne  current  model  is  very  similar  to  that  used  to  find 

the  direct  detection  current  model  of  this  thesis.  Lyons' 

results  are  presented  below. 

The  heterodyne  receiver  current  is  (Ref  16:31) 

( ka^ 

ih C t: ) = A3  j a(a+vt)h(a)exp[j0(a+vt)]exp[-j^|-]PDF(a)da  (85) 

where  A^  = anc*  a(°0>  h(a)»  6(a),  and  Ppp(a)  are  the 

o 

same  as  defined  for  the  direct  detection  current.  Due  to  the 
defining  equation  (Ref  Eq.  (4))  for  the  heterodyne  receiver 
current,  the  current  of  Eq.  (86)  is  centered  at  the  optical 
frequency  fQ-fjp.  Also,  this  current  is  complex  in  contrast 


to  the  direct  detection  current  which  is  real  (Ref  Eq.  (34)). 
The  expected  value  of  the  heterodyne  current  is  (Ref  16:31) 

. E[ih(t)J  = 0 (86) 

and  its  covariance  is  (Ref  16:36) 

C.  (t,t')  = Aj  f f i>DF(oO  PDF(a')a(a+vt)a(a^vt')h(a)h(a') 

(a'+vAt)  dada"  (87) 

where  P^(Aa)  is  given  by  Eq.  (47)  and  Aa=ct-a". 

The  heterodyne  current  is  centered  at  the  frequency  fQ-fjp. 
Thus  even  though  its  mean  is  equal  to  zero,  this  does  not  imply 
that  the  mean  of  the  amplitude  of  the  current  is  zero.  The 
current  amplitude,  denoted  A(t) , is  found  by  writing  the 
current  as  the  sum  of  its  real  and  imaginary  parts.  It  is 

i,  (t)  = I + j I. 
h^  1 r J l 

= A(t)  cos  Y(t)  + j A(t)  sinY(t)  (88) 

In  polar  form  this  becomes 

ih(t)  = A(t)  ejYCt)  (89) 

The  current  is  physically  broken  into  its  real  and  imaginary 

parts  by  the  well  known  quadrature  model  (Ref  24:238)  as  shown 

in  Fig.  8a.  Note  that  the  quadrature  outputs  are  the  real  and 

imaginary  parts  of  the  currents.  If  the  quadrature  outputs  are 

squared  and  then  added  together,  the  resulting  output  is  the 
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current  amplitude  squared,  A (t) , as  shown  in  Fig.  8b. 


Low  A(t)  cosy(t) 

Pass  — 

Filter 


A(t)  cos  ( 2irf  y p+y  ( t ) ) co s 2 it  £ j pt 


sin2ir  f t 
I F 


Low  A(t)  sinY(t) 

Pass  

Filter 


Fig.  8a  The  Quadrature  Model  16:40) 


A(t)cosY(t) 
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(•r 

A(t) sin  Y(t) 


Fig.  8b  Determination  of  A (t)  from  the 
Quadrature  Outputs  (Ref  16:43) 


Fig.  8 Determination  of  A (t) 
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The  mean  (Re f 16:43)  and  the  covariance  (Ref  16:45)  of 


A (t)  are 


and 


E [A  (t) ] = C.ft.t--  t) 


C 2 (t.O  = [C.(t,0] 
1 


(90) 


(91) 


By  substituting  Eq.  (87)  into  Eq.  (90),  the  mean  of  A^(t)  is 
expressed  in  terms  of  the  reflectance,  a(a)  , and  the  system 
function,  h(a)  as  follows: 

E [A2  (t ) ] = A2  J j h(a)h(a^)PDF(a)PDF(a')a(a+vt)a(a'+vt) 


P£ (Aa) dada ' 


(92) 


As  in  Eq.  (74),  P£(Aa)  can  be  approximated  by  the  Dirac  Delta 
function 


P£(Aa)  = A4  6 (Act) 

where  A^  is  the  area  of  P£(Aa)  which  is  given  by 
A4  = | Pj^  (Aa)  dAa 


(93) 


(94) 


Eq.  (92)  then  becomes 


E [A2 (t) ] = A^A4 


f h2(a)P2p(a)a2(vt-a)  da 


(95) 


which  by  a change  of  variables  is  written 

E [A2 (t) ] = A2A4  I a(x2)h2(vt-x)P2F(vt-x)dx 


(96) 
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Thus  the  mean  of  A (t)  can  be  represented  by  a linear  system 

2 2 2 
model  of  a (x)  convolved  with  h (x)P~p(x). 

2 

The  covariance  of  A (t)  is  found  by  substituting  Eq . (87) 
into  Eq.  (91) . It  is 


c ,•  (t,t') 

A^ 


h(vt-x)h(vt'-x)a(x)a(x') 


PJi(x-x')PDI.(vt'x)PDF(vt'-x')dxdx']  2 


(97) 


Thus  as  in  the  direct  detection  case,  the  double  integral  in 
the  above  equation  is  identical  to  the  output  correlation  of 
a filter,  h(x) , whose  input  is  a zero  mean  noise  process  with 
the  correlation  function 

^(x.x')  = a(x)a(x')P£(x-x')  (98) 

Once  again  the  square  of  the  integral  simply  means  that  the 
noise  can  be  represented  by  the  product  of  two  identically 
distributed,  statistically  independent  noise  processes.  The 
second  moment  model  for  the  amplituded  squared  of  the  heterodyne 
receiver  current  is  shown  in  Fig.  9.  As  with  the  direct 
detection  current  model,  this  model  can  easily  be  modified  to 
include  quantum  noise  effects.  These  noise  terms  are  included 
by  simply  adding  them  to  the  speckle  noise  term  described  by 
Eq.  (97). 

The  SNR  for  the  heterodyne  current  model  is  found  by 
substituting  Eqs.  (92)  and  (97)  into  Eq.  (82).  The  Resulting 
SNR  is 

SNR.  = 1 (99) 


Fig. 9 The  Model  for  the  Amplitude  Squared  of  the 
Heterodyne  Receiver  Current  (Ref  16:46) 


The  model  for  the  amplitude  squared  of  the  heterodyne 
current  and  the  direct  detection  current  model  are  very  similar. 
The  only  differences  are:  (1)  the  covariance  of  the  heterodyne 
current  model  has  a specific  value  rather  than  a bound,  (2) 
the  system  function  of  the  heterodyne  current  model  is  dependent 
on  both  h(x)  and  P^p(x) , and  (3)  the  two  models  differ  depending 
on  the  constants  Ap , A^ , A^ , and  A^ . Since  there  are  no 
major  differences  in  the  two  models,  the  comments  made  about 
the  direct  detection  model  are  valid  for  the  heterodyne  model. 
Thus,  as  the  heterodyne  system  function,  h'(x)  = h(x)PDp(x), 
becomes  wider  in  space,  the  higher  spatial  frequency  content 


of  a^(x)  is  lost  and  the  worse  the  resolution  ability  of  the 

2 

system  becomes.  The  model  of  E [A  (t)]  is  also  identical  to  the 
result  discussed  by  Goodman  (Ref  13:107)  for  an  incoherent 
imaging  system.  Thus,  square  law  detection  can  also  be  thought 
of  as  an  incoherent  imaging  system  with  an  additive  noise  term. 
Finally,  the  noise  in  the  heterodyne  model  is  also  signal 
dependent.  Thus  the  SNR  can  not  be  improved  by  increasing  the 
signal  power. 

One  important  result  in  the  comparison  the  two  models  is 
that  the  SNR  for  the  direct  detection  model  is  always  equal 
to  or  larger  than  the  SNR  of  the  heterodyne  model.  Tnerefore 
at  its  worst,  the  direct  detection  model  performs  as  well 
as  the  heterodyne  model.  Thus  on  the  basis  of  the  SNR,  the 
direct  detection  receiver  produces  the  better  results  for  the 
laser  line-scan  imaging  system.  However,  the  SNR  is  only  one 
of  many  parameters  used  to  evaluate  system  performance  and 
should  be  considered  as  only  a part  of  the  whole  picture. 


IV.  The  System  Models  for  Quasimonochromatic 
and  Broadband  Laser  Sources 


In  Chapter  III,  the  system  models  for  a monochromatic  laser 
were  determined.  These  models  are  now  extended  to  find  the 
resulting  models  for  broadband  and  quasimonochromatic  sources. 
First,  the  monochromatic  laser  field  is  extended  to  broader 
bandwidths.  Then,  by  using  the  new  laser  sources  and  following 
the  procedure  used  to  determine  the  previous  models,  the  current 
models  for  extended  laser  sources  are  found. 


quas  imonochromat ic  and  Broadband  Fields 

The  laser  field  given  in  Eq.  (23)  is  valid  for  monochromatic 
lasers  only.  This  field  can  be  modified,  so  that  it  is  valid 
for  both  quasimonochromatic  and  broadband  lasers,  by  a method 
known  as  linear  frequency  modulation  or  chirping.  The  theory 
of  chirping  is  discussed  in  a paper  by  Klauder,  et.  al. 

(Ref  25).  The  following  paragraph  is  a brief  summary  of  the 
results  of  chirping. 

Consider  a rectangular  time  pulse  with  width  T seconds  and 
denoted  rect[t/T].  The  bandwidth  of  this  pulse  is  Found , by 
Fourier  analysis,  to  be  approximately  yHz.  If  this  pulse  is 
chirped,  its  resulting  waveform,  denoted  UQ(t),  is  (Ref  25:754) 


UQ(t)  = rcct[t/T]  exp[jirAf  t ] 


(100) 


where  exp[jirAf  t ] is  the  chirp  waveform  and  Af,  whose  value 
is  chosen  by  the  system  design,  has  units  of  inverse  seconds 
squared.  It  is  shown  (Ref  25:  755-758)  by  Fourier  analysis 
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UQ(t)  has  a bniutwidt.li  of  ACT  - 1 Kllr. . Thus  t ho  bandwidth 

of  a signal  lias  boon  extended  by  an  order  of  magnitude  by 

(>  - 2 

chirping.  It  At  were  chosen  to  be  10  (sec)  , the  bandwidth 
would  have  increased  by  two  orders  of  magnitude.  The  point 
here  is  that  a wave  of  any  desired  bandwidth  can  be  produced 
by  chirping. 

This  point  also  holds  for  the  monochromatic  laser  field 
in  liq.  (23).  liq.  (23)  is  spatially  dependent  only.  Thus 
when  pictured  in  time,  it  has  a constant  amplitude  and  extends 
for  all  time.  Ideally  it  could  be  thought  of  as  a rectangular 
t i me  pu 1 se  w i t h i n f i n i t e width.  Howe v e r , i n reality  t he 
width  is  only  the  time  interval  for  which  the  laser  is  turned  on. 
In  the  laser  line-scan  imaging  system,  this  interval  could  be 
thought  of  as  the  time  needed  to  make  a single  scan  of  the 
object  surface.  In  any  case,  multiplying  liq.  (23)  by  the 
chirp  waveform  results  in  a field  which  is  not  monochromatic, 
but  rather  it  has  some  bandwidth.  This  bandwidth  i s determined 
by  Af  atul  the  t ime  interval  for  which  the  laser  is  on.  It 
is  reasonable  to  assume  that  the  chirped  laser  field  can  be 
designed  to  be  either  quns imonochromat ic  or  broadband.  In 
cither  case,  the  general  form  of  the  chirped  laser  field  is 


U(o.t) 


= — exp[j  (kz-J)  ]exp[j-^-a2] 

(Xz) 2 4 22 

2 7 

exp  [ - — ] exp[jmAf  t2]  (101) 

WZ(z) 

With  the  use  of  this  chirped  field,  it  is  possible  to  determine 
the  system  models  for  the  extended  bandwidth  sources. 

The  Quas imonochromat ic  System  Models 

A quasimonochroinatic  laser  source  is  used  with  both  a 
direct  detection  and  a heterodyne  receiver.  The  second  moment 
models  for  the  two  systems  are  determined.  First,  the  reflected 
laser  field  that  is  incident  on  the  detector  surface  is  found. 


The  Quas imonochromat ic  Field  Incident  on  the  Detector. 
Quas imonochromat ic  fields  are  defined  in  Appendix  A as  those 
which  suffered  a negligible  time  delay  during  propagating. 
Thus  quas imonochromat ic  fields  are  similar  to  monochromatic 
ones  in  that  they  both  propagate  independent  of  time.  For 
the  quas imonochromat ic  case  the  field  at  the  detector  is 
given  by  Hqs . (18)  and  (21): 


U,,(x,t)  = - 1-  i:  exp[-j  (kz-J)  ]exp[-jk|^] 


(Az) 


Ur.Qta.t) 


2 x« 

e*P[-jjY  a ] exp  [j  2*^-]  da 


(102) 


where 


U n(a,t)  = a(a+vt)exp[j0(a+vt)]  — — v exp[-  ] 

r’^  (Xz)*  r-(z) 

exp[j(kz-  J)  ]cxp[- j^a2]exp[-jTTAf  t2] 


(103) 
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Substituting  Eq.  (103)  into  Eq.  (102),  the  field  at  the 
detector  is 

2 

Ud(x,t)  = ^ exp[-j2(kz-  j)  ]exp[-jk|Y]exp[-jirAf  t“] 

| a (a+vt)  exp  [j9(a+vt)]h(a)exp[-j^-a^] 

exp [ j 2 *£]da  (104) 

Note  that  except  for  the  chirp  term,  Eq.  (104)  is  identical 
to  the  field  incident  on  the  detector  for  the  monochromatic 
case  (Ref  Eq.  (28)).  Eq.  (104)  is  now  used  to  determine  the 
direct  detection  model. 

The  Quasimonochromatic-Direct  Detection  Current  Model. 

The  ideal  direct  detection  receiver  current  is  given  by 
Eq.  (5)  as 

id  = $ T~  | j I Ud  (x  ,y)  I 2 dxdy  (5) 

In  Eq.  (5),  the  current  is  proportional  to  the  magnitude 

of  the  field  incident  on  the  detector  surface.  Taking  the 

magnitude  of  Eq.  (104)  cancells  the  chirp  term  with  its 

conjugate.  Thus  the  magnitude  of  the  field  incident  on  the 

detector  is  the  same  for  both  the  monochromatic  and  quas imonochromat i 

laser  sources.  The  result  is  that  the  second  moment  models 

for  the  direct  detection  current  are  identical  in  both  cases. 

This  finding  in  not  a complete  surprise.  From  Eq.  (20) 

and  for  a typical  value  of  I r „ , I = 10'^  meters,  the  quasi- 

monochromatic  field  must  have  a bandwidth  much  less  than 
5 

3x10  Hz.  TJius  with  respect  to  the  laser  frequency,  the 
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quasimonochromatic  field  is  almost  monochromatic  and  the  two 

direct  detection  current  models  should  be  similar.  In  this 

thesis  it  has  been  assumed  that  lasers  can  be  monochromatic, 

when  in  reality  those  lasers  referred  to  as  monochromatic 

actually  have  some  narrow  bandwidth.  The  above  results  shows 

that  assumption  of  modeling  a narrowband  source  as  monochromatic 

is  valid.  Finally,  since  both  direct  detection  current  models 

are  identical,  they  have  the  same  SNR  and  resolution  characteristics. 

The  heterodyne  current  model  is  now  investigated. 


The  Quasimonochromatic-Heterodyne  Current  Model.  The 
ideal  current  of  the  heterodyne  receiver  is  given  by  Eq.  (4) 
as 


ih(t)  = ElP  f Re{Ud(x,y)U*0(x,y)exp[- j2nf1Ft]  }dxdy  (4) 
o 


Substituting  Eq.  (104)  into  Eq.  (4),  the  current  is 


ih,Q('t')  = A3Re^exP  [~  j 2 (kz-  J)  ]exp  [- jirAf  t2]exp  [- j 2irfIpt] 

; | PD(x)a(a+vt)h(a)exp  [je  (a+vt)  ]exp[-j|a2] 

exPr-j£x2]UL*(x,t)exp[j2ir^]dadx}  (105) 

2 k 2 

In  the  above  equation,  the  exp  [-jirAf  t ] and  exp[-jyj  x ] terms 
can  be  cancelled  by  indentical  phase  terms  in  the  local 
oscillator.  The  x integral  is  now  the  spatial  Fourier  transform 
of  the  aperture  function  denoted 


Fx[pD(x)] 


PppC01) 


(106) 


Only  the  amplitude  of  the  current  is  of  interest,  thus  the 
phase  terms  in  front  of  the  integral  are  neglected.  Also, 

k 2 

exp[- j — ] = 1 as  shown  in  the  approximation  of  Eq.  (48). 

The  current  is  now  written  as 

ih  = A3  | PDF(a)a(a+vt)h(a)exp[j0(a+vt)]da  (107) 

This  result  is  identical  to  the  complex  current  of  Eq.  (85)  . 

Thus  as  resulted  for  direct  detection  current  models,  the 
heterodyne  current  models  are  identical  for  the  monochromatic 
and  quasimonochromat ic  laser  sources.  The  comments  made  at 
the  end  of  the  previous  section  for  the  direct  detection 
current  models  could  also  be  stated  for  the  heterodyne  current 
models . 

The  Broadband  System  Models 

The  second  moment  models  are  found  for  two  systems:  a 
broadband  laser-heterodyne  receiver  and  a broadband  laser- 
direct  detection  receiver  system.  The  laser  field  is  assumed 
to  be  broadband  due  to  an  appropriate  chirp  term.  First,  the 
field  incident  on  the  detector  is  determined. 

The  Broadband  Field  Incident  on  the  Detector.  A broadband 
field  was  defined  in  Appendix  A as  one  which  suffers  a time 
delay  during  propagation.  Because  of  this  delay,  the  mathematics 
involved  in  propagating  a broadband  field  are  somewhat  more 
difficult  than  for  a monochromatic  or  quasimonochromatic  field. 

The  configuration  for  determining  the  broadband  field  is  illustrated 


10  Notation  for  the  Broadband  System 


Consider  the  setup  of  Fig.  10.  The  monochromatic  field 


at  the  laser  output  is  given  by  Siegman  as  (Ref  15:307) 


where  Wq  is  the  spot  size  of  the  laser  beam  at  the  laser  output 
The  broadband  laser  field  at  the  laser  output,  denoted  U.  R(B) 


The  broadband  field  incident  on  the  rough  surface,  denoted 


U-  D ( a)  , is  found  by  propagating  the  laser  field  by  the  extended 
i »B 

Huygens -Fresnel  integral  given  in  Eq.  (18)  . The  incident 
field  is 


Ui,B(o°  = 


(Az) 


Tj  exp[j  (kz-J)  ]exp[jk^-] 


exp [ - ] 

W. 


exp  [ j tt  A f (t |^)2]exp[j||  ]exp[-j^|^-  ] da  (110) 


4 -i 

where  01  is  the  propagation  time  delay  from  the  laser 
c 

to  the  surface.  The  field  reflected  from  a rough  surface 
was  given  in  Eq.  (25).  The  reflected  broadband  field  is 


Ur  R(a,t)  = a(a+vt)exp[j9(a+vt) ]Ui>B  (a,t) 


(111) 


The  field  incident  on  the  detector  is  determined  once  again 
using  the  extended  Huygens  - Fresnel  integral.  It  is 


1 v 2 

Ud,B(x’t}  = 77-1?  lexP[i  (kz-J)]exp[-jg-] 


(AZ)  1 


ri~J  j 


Ur(a,t- 


exP  t * ] exp  d 


(112) 


where  r is  the  progation  time  delay  from  the  surface  to 
the  detector.  Substituting  Eqs.  (110)  and  (111)  into  Eq.  (112) 


A v x2 

^ exp  [j(kz-^-)]exp[-j  ky^]  a (a+vt)  exp  [ j 0 (a+vt)  ] 


exP  [' j-^— ] exp  [-—  2]exp(- j^-  ]exp[j  j|(aB+xa)  ] 


exp[-jirAf(t  - 


-2i  )2]  dadB 


CH3) 


From  Fig.  2 and  the  approximation  of  Eq.  (11)  it  is  seen  that 


rn1  = [z  + Cot-B) 


1 * TT  (“-S)2 


(114) 


[z2  + (x-a)2] 


z * (x-a)2 


(115) 


By  use  of  the  last  two  equations,  it  can  be  shown  that  the  exponent 
in  the  chirp  term  is 


r r 

j*Af(t-  -|i  - )2  = jirA£[d1(t)  + d2(t){(a-g)2  + (x-a)2} 


+ - 2L2-  (a-B)2(x-a)2] 

2z  c 


(116) 


where 


dx(t)  = 


t2  + 


(117) 


d2^  -jr-ir 


(118) 


The  last  term  in  Eq.  (116)  corresponds  to  the  third  term  in  the 
binomial  expansion  given  in  Eq.  (11).  Thus  it  can  be  neglected. 
After  combining  Eq.  (116)  with  Eq.  (113)  , the  field  at  the 
detector  is 


Ud  B(x,t)  = ^ exp[-j2(kz-J)]exp[-jk^-]exp[-jir&fd1(t)  ] 

2 2 

a(a+vt)  exp  [j  0 (a+vt)  ] exp  [ - j^|—  ]exp(-j^-  ] 


exp  [*— 2 1 exp  [j^ytag  + xa)  ] 
o 


exp  [- jAfird2  (t)  { (a-$)  ^ + (x-a)2}]  dadB 


(119) 


The  double  integral  in  the  a bove  equation,  denoted  Q, 
can  be  written 


Q = ||  g1(a,g)exp[jAf  g2(a,0)]dadB 


(120) 


where 


gjfa.B)  = 


a (a+vt) exp [ j 9 (a+vt) ]exp [-j— — ] 


exp  [- j^-]  exp  [-— 2 ] exp  [j~(a0+xa)  ] 
Wo 


(121) 


g2(a,e)-  = “ tt  d 2 (t)  [ (a-B)  2-  + (x-a)2] 


(122) 


- l \ 


Eq.  (120)  has  the  form  of  an  integral  which  can  be  solved  by 
the  method  of  stationary  phase.  This  method  is  discussed 
in  Appendix  B and  can  be  used  to  show  that  as  the  quantity  Af 
becomes  large,  Q can  be  approximated  by 

2 

Q = ( ~ j A f d 2 (t ) ] ~ 1 a(x+vt)exp[j0(x+vt)  ]exp[j^-] 


exp [-  — ] 

"o2 


(123) 


It  is  emphasized  that  Eq.  (123)  is  only  an  approximation  which 
becomes  more  exact  as  Af  becomes  larger.  No  attempt  is  made 
in  this  paper  to  determine  the  relationship  between  Af  and 
the  variation  of  Q from  its  actual  value.  However,  this 
relationship  has  been  analyzed  by  Cook  and  Bernfeld  (Ref  27: 
Chapter  3)  and  by  Fowle  (Ref  28:  61-67)  for  integrals  similar 
to  Eq.  (120) . Thus  it  seems  feasible  that  Eq.  (123)  can  be 
quantified  by  an  analysis  similar  to  the  methods  of  above  two 
references . 

Substituting  Eq.  (123)  into  Eq.  (119)  yields  the  field 
incident  on  the  detector  which  is 


Ud,B^x,t^  = Af  ( 2z~  tcY  exp[-j2kz]exp[-j7rAfd1(t)] 

x2 

a(x+vt)  exp [ j 0 (x+vt) ] exp [- — 7 ] 

W 

o 


(124) 


This  equation  is  now  used  to  determine  the  receiver  currents. 

The  Broadband  Direct  Detection  Current  Model.  The  ideal 
direct  detection  receiver  current  is  given  by  Eq.  (5).  Combining 
Eq.  (5)  with  Eq.  (124)  produces  a broadband  current  of 
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(125) 


1d,B 


(t)  " A, 


a (*)  h (.vt  - x)  dx 

U j . 


where 


A 


4 


AXfo^  2 

I-Af  (2z-tc)  ^ 


(126) 


and 

2x^ 

h (x)  = exp  [ -— j]  Pn(x) 
0 VI_  u 


(127) 


hQ(x),  like  h(x) , is  also  called  the  system  function. 

Note  that  the  broadband  current  of  Eq.  (125)  contains 
no  random  terms.  This  is  an  amazing  fact  because  it  implies 
that  there  is  no  speckle  noise  present  in  the  current.  Most  laser 
line-scan  imaging  systems  are  limited  by  speckle  noise.  In 
fact,  the  four  systems  previously  investigated  in  this  paper 
had  signal  dependent  noise.  This  is  obviously  not  the  case 
for  this  system  since  it  has  no  speckle  noise.  Thus  from  the 
above  results,  it  appears  that  a laser  line-scan  imaging  system 
which  is  not  speckle  noise  limited  can  be  designed. 

2 

Eq.  (125)  has  the  form  of  a convolution  process  of  a (x) 
with  hQ(x).  Fig.  11  is  a diagram  of  this  model.  The  quantum 
noise  effects,  which  were  neglected,  can  be  included  by  adding 
the  appropriate  noise  terms  to  the  broadband  current,  i^  g(t). 

The  new  system  function,  hQ(x)  , is  dependent  on  the  size  of 
the  detector  surface  and  the  beam  spot  size  at  the  laser  output, 
WQ.  It  serves  the  same  purpose  as  the  previous  system  function, 


56 


Fig.  11  Broadband-Direct  Detection  Current  Model 


h(x).  It  was  shown  in  Eq.  (84)  that  as  h(x)  becomes  wider, 
the  system's  resolution  ability  becomes  worse.  This  comment 
is  true  for  hQ(x)  also. 

The  Broadband-Heterodyne  Current  Model.  The  broadband 
heterodyne  current  is  found  by  substituting  Eq.  (124)  into 
Eq.  (4).  The  latter  equation  gives  the  heterodyne  receiver 
output  current  which  becomes 


ih,B(t) 


= A,.  Re(exp  [- j2kz]exp  [- jirAfd^(t)  ] 

exp [- j 2nf j pt ] a (x+vt) exp [ j 0 (x+vt ) ] 

2 * 

PD(x)exp  [-— 2]Ulq  (x,t)dx 


(128) 


where 


a,  ■ Sa 


(129) 


In  Eq.  (128),  the  exp[-jTrAf  d^t)]  can  be  cancelled  by  an  identical 
term  in  the  local  oscillator  field.  Only  the  current  amplitude 
is  of  interest,  thus  the  phase  terms  in  front  of  the  integral 
can  be  neglected.  The  complex  current  is  now  written  as 


ih,B(t:i  ” A5 


PD(x)  a(x+vt)  exp [ j 9 (x+vt) ] 


exp  [-— 7 ] dx 


C130) 


W 


The  rough  surface  phase  term  is  the  only  random  term  in 
Eq.  (130).  Thus  the  mean  of  the  current  is  given  by 

°92  f 

E [ih,B^t^  = A5  exP[ T~2^  JPDW 

2 

x 

a(x+vt)  exp[-^“ 2]dx  (131) 

o 

2 

where  exp[-  ] is  the  value  of  the  characteristic  function 

from  Eq.  (41).  For  a surface  with  oh  _>  X in  Eq.  (28),  the 
value  of  the  characteristic  function  is  approximately  zero 
(Ref.  16:31).  Thus  the  mean  of  the  current  is  effectively 
equal  to  zero; 

E[ih,B(t)1  = ° (132) 


The  covariance  of  the  current  is 
Elih,B(t)  1h,B(t")1  = A5 


PD(x)  PD(x')a(x+vt) 


a(x'+vt')  hj(x)  hj(x') 
P^(Ax+vAt)  dxdx ' 


(133) 


where  h. (x)  = exp[-  — ?]  • Because  of  the  zero  mean,  the 
1 

correlation  of  the  current  is  equal  to  the  covariance. 

Comparing  Eqs.  (131)  and  (132)  to  Eqs.  (86)  and  (87)  reveals  that 
the  broadband  heterodyne  current  model  is  almost  identical  to 
the  monochromatic  model.  The  only  differences  in  the  two  are  their 
amplitudes  and  their  system  functions.  These  are  minor  differences. 

Both  system  functions  serve  the  same  purpose;  they  determine 

the  system  resolution  ability.  The  broadband  system  function, 

h^(x),  has  a width  dependent  on  the  beam  size  of  the  laser  output, 

but  the  monochromatic  system  function,  h(x) , has  a width  dependent 

on  the  beam  size  at  some  distance  z from  the  laser.  Because  the 

beam  size  widens  as  the  field  propagates  further  from  the  laser, 

h(x)  is  wider  in  space  than  h^(x).  Thus  the  broadband  heterodyne 

system  will  have  better  resolution  characteristics  than  either 

the  monochromatic  or  quasimonochromatic  heterodyne  systems. 

All  three  systems  have  the  same  SNR's.  Comparing  h^(x)  to 

hQ(x),  the  broadband  direct  detection  system  function,  shows 
2 

that  hQ(x)  = h^(x).  It  can  be  shown  by  Fourier  analysis  that 
hQ(x)  is  n'rrower  in  space  than  h^(x).  Thus  the  broadband- 
direct  detection  system  has  the  best  resolution  ability  of 
all  the  systems. 

In  this  chapter,  four  system  models  were  found.  The 
quasimonochromatic-direct  detection  and  quasimonochromatic- 
heterodyne  current  models  were  shown  to  be  identical  to  their 
respective  monochromatic  models.  The  broadband  systems  were 
found  to  have  better  resolution  ability  than  the  narrowband 
systems.  The  major  discovery,  however,  is  that  the  broadband- 
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direct  detection  system  does  not  s iffer  from  speckle  effects 
and  it  has  the  best  resolution  ability  of  all  six  systems 
investigated. 
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The  Application  of  the  Laser  Line-Scan 
Imaging  System  to  the  Ranging  Problem 

In  this  thesis,  ranging  is  considered  to  be  the  process 
of  determining  the  height  of  some  point  on  the  rough  surface 
relative  to  other  points  on  the  surface.  This  process  measures 
macroscopic  variations  in  the  surface  (on  the  order  of  centimeters) 
as  opposed  to  the  microscopic  variations  in  the  surface  (on 
the  order  of  the  optical  wavelength).  The  laser  line-scan 
imaging  system  is  used  to  perform  the  process  of  ranging. 

In  this  chapter,  the  method  and  system  used  in  ranging 
is  discussed  first.  The  selection  of  the  intermediate  and 
carrier  frequencies  in  design  considerations  is  presented. 

Then  the  second  moment  model  due  to  speckle  effects  is  found 
and  quantum  noise  effects  are  added  to  the  model.  This  model 
is  then  used  to  evaluate  the  system  performance  for  ranging. 


where  f is  the  modulation  frequency.  During  the  propagation 
from  the  laser  to  the  rough  surface  and  back,  the  modulation 

r01 

field  suffers  a time  delay  of  2— . r ^ is  shown  in  Fig.  2 and 
given  by  Eq.  (12).  The  laser  can  be  designed  with  a bandwidth 
that  is  narrow  enough  so  that  r^^  = z.  This  means  that  any 
changes  in  the  laser  field  away  from  z are  small  enough  to 
be  ignored.  The  bandwidth  requirement  is 


1 


> > 


I a~x I max 
2zc 


(135) 


2 

The  time  delay  is  now  2 — . The  modulation  field  that  is 
incident  on  the  detector  surface  is 

md(t)  = [1  + cos  2TTfc(t  - ^ )]H  (136) 

This  result  is  equivalent  to  the  quasimonochromatic  case  of 
the  previous  chapter  in  that  the  modulation  term  is  independent 
of  spatial  variations. 

If  the  time  delay  term  of  Eq.  (136)  can  be  isolated, 
then  the  distance  from  individual  points  on  the  rough  surface 
to  the  scanning  system  can  be  determine.  The  system  used 
in  determining  this  distance  is  a quasimonochromatic  laser- 
heterodyne  receiver  system.  The  quasimonochromatic  laser 
is  chosen  because  like  an  actual  laser  it  has  a bandwidth, 
but  it  can  be  modeled  by  the  monochromatic  results.  The 
ranging  model  for  this  system  is  developed  in  the  following 
section. 

The  Ranging  Model 

As  with  the  other  models,  the  field  incident  on  the 
detector  is  found  first.  The  mean  and  covariance  of  the  current 
are  then  found.  Quantum  noise  effects  are  added  to  the  current. 
The  resulting  current  is  processed,  as  shown  in  Fig.  8,  to 
determine  its  amplitude  squared.  A process  is  then  developed 
to  isolate  the  time  delay  term  and  thus  allow  the  system  to 
be  evaluated. 


The  Modulated  Quasimonochromatic  Field  Incident  on  the 
Detector . The  modulation  term  is  dependent  on  time  and  the 


propagation  distance  z.  Thus,  this  term  is  not  effected  by  the 
Huygens-Fresnel  integral  equation  which  integrates  over  spatial 
coordinates.  It  can  be  thought  of  as  a constant  during  the 
propagation  process.  Therefore,  the  field  incident  on  the 
detector  is  the  same  as  for  the  monochromatic  case  (Ref  Eq.  (31)) 
but  now  it  is  multiplied  by  the  time  delayed  modulation  term 
of  Eq.  (136).  The  incident  modulated  field,  denoted  m(x,t) 
is 

Ud,mCx»t)  = mdCt}  UdCx) 

= [1  + cos  2-nrf c (t -~)  -A-exp  [- j 2 (kz-  J)  ] 

2 r 2 

exp[-jki-]  a (a+vt) exp  [ j 9 (a+vt)  ] exp  [ — ^ ] 

J tr(z) 

v 2 

exP  ["  j-^r— ] exp  [j  2iTyY  ]da  (137) 

By  algebraic  manipulation  it  can  be  shown  that 

[1  + cos27rfc(t-  ^ )]h  = (2)h  cos[2tt^  (t--^)]  (138) 

Therefore,  the  complex  envelope  of  the  field  in  Eq.  (137)  is 
centered  about  the  frequency  fc/2  (recall  that  the  optical 
frequency,  f , is  suppressed)  and  has  bandwidth  B„  which  is 
restricted  by  Eq.  (135).  A representation  of  the  temporal 
frequency  spectrum  of  Eq.  (137)  is  shown  in  Fig.  12a. 

The  mechanics  of  the  heterodyne  system  add  a local  oscillator 
field  to  the  incident  field.  The  complex  envelope  of  the  local 
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Ft[md(t)U,Cx)] 


Fig.  12a  Temporal  Frequency  Spectrum  of  the  Detector 
Field 


Ft[cos  2itf  tULQCx)] 


Fig.  12b  Temporal  Frequency  Spectrum  of  the  Local 
Oscillator  Field 


f f 

IF'-r 


f + ^c 
IF  — 
lt  2 


Fig.  12c  Temporal  Frequency  Spectrum  of  |U,  (x,t)  + 

UloCx.UI2 


Fig.  12  Temporal  Frequency  Spectrums  for  the  Heterouyne 
Receiver  Fields 


I 


oscillator  field,  denoted  U^fa)  > is  centered  at  the  intermediate 
frequency,  fjpi  aad  it  can  be  thought  of  as  having  some  narrow 
bandwidth,  Bq.  As  shown  in  Chapter  II,  the  heterodyne  current 
is  proportional  to  the  magnitude  squared  of  the  sum  of  the 
two  fields  which  is 

lmd(t)Ud^  + cos  [2TTfIFt]UL0(x)  I 2 = m^(t)  | Ud(x)  | Z 

+ (cos  2 TTf  J pt ) 2 |UL0(x)  I 2 

+ md(t)  cos(2irfIFt)2Re{Ud(x)ULO(x)  ) 

= 2 [1  + cos(27rfc(t  - ^)>]  | Ud(x)  |2 

+ [1  + cos(27r2fIFt)]  |UL0(x)  !2 

f f 

+ (2)p2[COs{2TT(fIF  + -±)  (t  - ^0)  + COS  { 2 Tr  f£j.  F - -±)  (t-^-)}] 

Re(Ud(x)U*0(x)}  (139) 

Fig.  12c  is  a representation  of  the  one  sided  temporal 

frequency  spectrum  of  Eq.  (139) . Recall  that  the  heterodyne 

system  has  a bandpass  filter  centered  near  the  frequency  fjp. 

Thus  the  heterodyne  current  is  dependent  on  only  the  last  term  in 

Eq.  (139).  However,  the  last  term  contains  two  signals;  one  is 

f _ ' 

centered  about  the  frequency  fTC  - — and  the  other  is  centered 
fc  IF  2 

about  fjF  + -j-  . Only  one  signal  is  needed.  Thus  it  is  arbitrarily 

I! 


MI* 
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chosen  that  the  signal  centered  about  fjp-  -y  is  filtered 

out.  This  can  be  done  with  the  sue  of  a bandpass  filter 

f 

with  a bandwidth  of  Bn  + B~  centered  at  fT„  + The 


remaining  signal  is  used  to  determine  the  heterodyne  current. 

The  signals  illustrated  in  Fig.  12c  may  overlap  with 

each  other  depending  on  the  values  of  f c , fjp»  Bp,  and  Bq. 

This  overlapping  is  called  aliasing.  To  prevent  aliasing 

f 

between  the  signal  centered  at  fjp  + -y  and  the  other  signals 
the  following  three  conditions  must  be  met: 


These  three  conditions  reduce  to  two  restrictions 


Fig.  13  is  a graph  of  the  possible  values  of  f and  fTC 

C £i  r 

which  avoid  aliasing  between  the  signal  about  fjp  + -y  and 
the  other  signals.  The  frequency  locations  in  Fig.  12c 
and  the  restrictions  of  Fig.  13  have  assumed  that  fTC  < f . 


Very  similar  results  would  have  been  obtained  if  it  were 


assumed  that  f > f 


values  which 
prevent  aliasing 


(£  +3Bn+Bn) 


Fig.  13.  Values  of  fTP  and  f which  Prevent  Aliasii 


For  the  restrictions  of  Eqs.  (143)  and  (144)  and  for  the 

f 

bandpass  filter  centered  at  fJF  + , the  resulting  field  is 


UBp(x,t)  = (2)  cos[27r(fIF  + -^)(t-^)]Re(Ud(x)UL0(x)}  (1 


This  field  is  centered  about  fjF  + and  has  bandwidth  BR  + Bc 

In  this  case  the  value  of  intermediate  frequency  is  unimportant 

Thus  by  appropriate  modulation  and  filtering  techniques,  the 

f 

center  frequency  of  the  field  can  be  changed  to  . 


6 


The  resulting  field  still  has  bandwidth  BR  + BQ,  and  is 

UBp  (x,t)  = (2)**  cos[2ir-^  (t-  ^-)]Re(Ud(x)U*0(x)} 

= Re{Ud>m(x,t)U*0(x)}  (146) 

where  Ud  m(t)  was  given . in  Eq.  (137).  This  field  is 
now  used  to  determine  the  heterodyne  receiver  current. 

The  Ranging  Current  Model.  The  heterodyne  current  is  found 
by  substituting  the  bandpassed  field  of  Eq.  (146)  into  Eq.  (4) . 
The  result  is 


HF1  TzT  f1  + cos27rfc(t- 


Re( 


exp [- j 2 (kz-  ?)] 


2 2 

exp(-jky-]a(a+vt)exp[— f ] 

^ VT(z) 


PD  (x)  exp  [ - j k^-]  exp  [ + j ULQ  (x)  dxda 


(147) 


As  with  previous  heterodyne  models,  the  exp[-jk^-]  term  can  be 
cancelled  by  an  identical  phase  term  in  the  local  oscillator 
field.  Then  the  x integral  is  just  the  Fourier  transform  of 
the  aperture  function  denoted 


W*)i 


fx  • 


a 

Xz 


= PDF(a) 


(148) 


Also,  only  the  amplitude  of  the  current  is  of  importance  here, 
so  the  exp [ - j 2 (kz - j)]  phase  term  can  be  neglected.  The  current 
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Eq.  (5)  with  Eq.  (124)  produces  a broadband  current  of 


is  then  written  as 


W1*  ■ V1  * cos2„fc(t-f)]'S  [ Pjjp(a) 


a(a  + vt)h(a) exp [j  0 (a+vt ) J exp [- j — ] da 


(149) 


Comparing  Eq.  (149)  to  Eq.  (85),  it  is  seen  that  the  modulation 
heterodyne  current,  except  for  the  time  delayed  modulation  term, 
is  identical  to  the  monochromatic  case  and  thus  the  quas imonochromat ic  I 
case.  This  term  is  not  random  and  its  amplitude  varies  between 

j, 

zero  and  (2)  2 , thus  it  does  not  effect  the  mean  and  covariance 
of  the  current.  They  are  the  same  as  before: 


^h.m^l  = 0 


(150) 


C.  (t) 
1h,m 


md(t)md(t')  C.  (t,t-) 
h 


2 2.1  n ^ 


2 Z s ■,  •,  *5 


A3  [1  + COs(2TTfc(t-^)}]  2[1  +COS{2TTfc(t'-  } ] 


PDF PDF a (a  + vt ) a(a'+vt>)  (Aa+vAt) 


h(a)h(a") dada^ 


(151) 


To  make  the  ranging  problem  more  general,  quantum  noise 
effects  are  now  included.  Quantum  noise  is  denoted  n-p^t).  It 
is  included  by  simply  adding  it  to  the  modulated  current  of 
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Eq.  (149).  This  sum,  denoted  iT(t)  for  total  current,  is 

iTCt)  = md(t)ih(t)  + nTH(t)  (152) 

The  first  two  moments  of  the  quantum  noise  due  to  heterodyne 
quantum  effects  are  (Ref  26:176) 


E [nTH(t)  ] = 0 

(153) 

<l2nA, 

cn  (t.t-)  = a-  6(t-0 

TH  hf 

(154) 

Therefore,  the  mean  and  covariance  of  the  total  current  are 
E [iT(t) ] = E (ih>m(t) ] + E[nTH(t)j 

= 0 (155) 

and 

C.  (t,t')  = m,(t)m,(t")C.  (t,t')  + C (t,t')  (156) 

h 1h  nTH 

It  can  be  seen  from  Eq.  (154) , that  quantum  noise  is  white 

noise.  That  is,  its  power  spectral  density  has  a constant 

amplitude  and  exists  for  all  frequencies.  One  way  to  decrease 

the  power  in  the  noise  is  to  pass  the  noise  through  a bandpass 

filter.  The  total  current  can  be  passed  through  a bandpass  filter 

f 

with  bandwidth  BR  + BQ  centered  about  — . This  does  not 

effect  the  modulated  heterodyne  current  but  it  does  decrease  the 

quantum  noise.  The  bandpassed  quantum  noise  is  denoted  n„„  _n(t) . 

Trl . BP 
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If  f > > Bd  + Bn  (possible  by  system  design) , the  quantum 
C K U 

noise  can  be  represented  in  terms  of  its  quadrature  components 
(Ref  24:237).  The  quantum  noise  becomes 


nTH  gp  (t)  = nc(t)  cos  (2tt—  t)  -ng  (t)s  in(27r-^  t) 


= r (t)  cos(2ir--£  t + <j>TH(t)) 


(157) 


where 


r(t)  = [n^ (t)  +n^(t)]2 


(158) 


U)  ^ 


-1  rns(t) 


(159) 


The  quantum  noise  still  has  zero  mean. 

By  definition  the  covariance  of  the  bandpassed  quantum 
noise,  now  denoted  CTH(t,t'),  is 


cTH(t»t') 


E ^nTH,BP('t-*  NTH,BP('t^  ^ 


f f 

E [cos  { 2ir-^-t  + ^p^ft)  }cos(2ir-y-t"  + i$>TH(t')} 


r (t) r (t ')  ] 


(160) 


The  covariance  of  the  total  current  including  bandpassed  quantum 


noise  is 


Ci  (t,t')  = md(t)md(t')Ci  (t,t')  + CTH(t,t') 
T h 


(161) 
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Fig.  14  Model  of  the  Total  Current 

The  noise  terms  represented  by  the  total  current  covariance 
are  shown  in  Fig.  14. 

In  Eq.  (163),  the  covariance  term,  C.  (t,t'),  is  not 

xh 

stationary  (Ref  Eq.  (98)).  For  simplification  it  is  assumed 

that  the  surface  reflection  coefficient,  a(x)  , is  a constant. 

This  implies  that  C.  (t,t^)  is  stationary  and  thus  Fourier 

xh 

techniques  can  be  used  to  examine  this  term.  Note  that  since 
the  quantum  noise  is  stationary,  then  C^,^(t,t')  is  stationary 
also.  Eqs . (162)  and  (163)  are  used  to  determine  the  second 
moment  model  for  ranging. 

The  amplitude  squared  of  the  total  current  is  found  by  the 
process  shown  in  Fig.  8.  The  first  two  moments  of  the  amplitude 
squared  are  given  by  Eqs.  (90)  and  (91).  They  are 

E(A^(t)]  = C.  (t,t--  t)  • 


systems.  The  major  discovery,  however,  is  that  the  broadband- 


Eq.  (162)  is  the  "signal"  component  of  the  model.  The 

first  term  in  this  equation  is  composed  of  two  signals: 

one  is  centered  about  zero  frequency  and  the  other  about 

2 

The  amplitude  of  these  tvro  signals  is  C.  (t,t'  =t)  or  a.  (t) 

1h  xh 

The  last  term  in  Eq.  (162)  is  the  variance  of  the  bandpassed 
quantum  noise.  This  term  is  centered  about  zero  frequency 
also  (Ref  24:239).  If  the  "signal"  terms  are  A.C.  coupled, 
the  terms  at  zero- frequency  are  filtered  out.  This  leaves 


Eq.  (163)  describes  the  "noise"  components  of  the  model 
Let  the  first  term  in  this  equation  be  denoted  n,-.(t).  It 


describes  a noise  process,  call  it  n (t) , which  is  zero  mean 


i 


has  a covariance  of  [C.  (t,t 'll",  and  is  multiplied  by 

_ * h 

[1  + cos(2irf.(t  - ~) ) | . Multiplication  of  n (t)  by  1 yields 

C **  * I 

a zero  frequency  term  which  can  also  be  filtered  by  A.C. 
coupling.  fig.  ISa  is  a representation  of  n^Ct). 

The  second  term  in  bq.  (163")  describes  two  statistically 
independent  noise  processes,  n.7(t)  and  n j - f 1 1 , which  arc 
multiplied  together.  nj7(t)  has  zero  mean,  a covariance  of 

2 T 1. 

C.  (t,t  ),  and  is  multiplied  by  fl  + cos { 2nf  (t --—)}]* . n,,(t) 

'li  c c 

lias  zero  mean,  a covariance  of  C.  (t,t'),  and  can  be  represented 

1 h r 

by  the  multiplication  of  r(t)  by  cos(2n-v  t + <K||(t))  (.Kef  I:.q.  (1 
Thus  the  multiplication  of  i\j,(t)  and  n ^ 1 1 ) , denoted  n^ft), 


n j j ( 1 ) * n j 2 ( t ) n ^ ( t ) 

■>„  L f 

- [1  + cos(2iri\,(t  *“)  1 l^cos(2it-y  + *TH(t)  ^ (t)  r (t  1 

* cos[2nfc(t  -“•)  + tTn  f 1 ^ 1 n l 5 1 r ^ 

+ c.os[2irfc  £ + <fTH(t)  liij  5Ct)  (1651 

where  n,,(.t)  * ( Vd  ' n,,(t)r(t).  The  last  term  in  bq . (165)  is 
at  zero  frequency  since  cos(2Trf  -•  * ^.^(t))  is  just  some  phase 
shift.  Thus  as  before,  ibis  term  can  be  filtered  out.  The 
remaining  term  is  represented  in  Tig.  15b. 

The  last  term  in  Tq . (163),  can  be  represented  by  two 
statistically  independent  noise  processes,  n^it)  and  tij7(t). 
bach  of  these  is  zero  mean,  has  a covariance  of  (. t , t *)  , and 
can  be  represented  by  the  product  of  r(t)  and  cos(2n  , t + ^...(tl 


Fig.  15a  Speckle  Noise  Representation 


Fig.  15b  Speckle-Quantum  Noise  Representation 


Fig.  15c  Quantum  Noise  Representation 


The  product  of  n^(t)  and  n^(t)  is 


n18(t) 


nl6(t)nl7(t) 

f f 

r1Ct)r2(t)cos(2iT-—  t + 4^(0)  cos(2tt-^  t + <J>7(t)) 
[cos  { 27rf  ct  + 4>1  (t)  + <Pz(t)} 

+ cosU^t)  - 4>2(t)H  aTH(t)  (166) 


where  n^tt)  = ^-r^  (t)  r9  (t)  . Once  again  the  last  term  in  Eq.  (166) 
is  centered  about  zero  frequency  since  cos(<j>^(t)  - 4>?(t))  is 
only  a phase  shift.  Thus  this  term  is  filtered  out.  The 
remaining  term  is  shown  in  Fig.  15c. 

A representation  of  the  "signal"  and  "noise"  erms  described 
by  Eqs.  (162)  and  (163)  has  now  been  suggested.  Tnus  from 
a second  moment  model  standpoint,  it  can  be  said  that  the 
araplituded  squared  of  the  total  A.C.  coupled  current  is 

A^AC(t)  - cos2„fc(t-  O^Ct) 

* cos2»fc(t  - f)  »sp(t)  | 

+ cos[27rfc(t  - ^)  + <J>TH(t)]  n]5(t) 

+ cos[2iTfct  + 4>  x ( t ) + <t>  2 ( t ) ] nTH  (t)  (167) 

Eq.  (167)  is  now  referred  to  as  the  ranging  current.  Fig.  16 
is  a diagram  of  the  second  moment  model  of  the  ranging  current. 


L 
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cos2tt£  (t 


cos2?rf  (t 


cos  [ 2tt f _ (t-— ) 


cos  [ 2tt f t+  <{>.  (t)  + <J>_  (t)J 


Fig.  16  The  Current  Model  for  Ranging 


With  the  help  of  Eq.  (167)  and  the  current  model,  it  is  possible 

2 z 

to  detect  the  time  delay,  — . This  is  done  in  the  following 
section. 


Detection  of  the  Time  Delay.  It  is  desired  to  measure  the 
time  delay  so  that  the  distance  from  the  rough  surface  to  the 
imageing  system  can  be  obtained.  What  is  actually  measured 
though  is  the  phase  delay  that  results  in  the  field  as  it 
propagates  down  and  back  up.  Measuring  the  phase  delay,  which 


2 7 , 

is  2 it f c — , is  good  enough  since  and  c are  known.  Thus, 
z,  the  distance  that  is  desired  to  be  known,  is  the  only 
unknown  and  can  be  easily  determined  from  the  phase  delay. 

However,  because  quantum  and  speckle  noise  are  present, 
the  phase  delay  that  is  measured  includes  the  phase  delay  of 
the  noises  as  well  as  the  laser  field's  phase  delay.  Therefore 
there  in  an  error  in  the  value  of  z that  is  measured.  The 
phase  delay  error  can  be  seen  from  a phase  diagram  of  the 
ranging  current  in  Eq.  (167).  First,  this  current  must  be 
rearranged. 

Consider  the  third  term  in  Eq.  (167).  It  can  be  rewritten 
as 

cos  [2rrfc(t-|)  + <f>TH(t)  ]n15  (t)  = cos  [2irfc (t-^r)  + ( 2tt f c | 

+ 4>TH(t))]n15(t) 

= cos2irfc(t-^-)cos  [2irf c |+<J>TH(t)  ] n15  (t) 

+ sin[27rfc(t  - ^-)  ]sin[27Tfc|+  4>TH ( t ) ]n15  (t)  (168) 

The  last  term  in  Eq.  (167),  can  be  rewritten 
cos  [ 2 tt f c t + <J>1  (t)  + d>2  (t)  ]n'H(t) 

= cos[2Trfc(t-^)  + 2Trfc^-  + ^(t)]  nJH(t) 

= cos  [ 2 tt f c ( t ■ ^-)]cos[2Trfc^-+  *3  (t)  ] n'H (t) 

+ sin[2irfc(t-^)]sin[2TTfc^  + (frjd)]  n{H(t)  (169) 
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where  4>3(t)  = 2 C t ) + Now  the  ranging  current  becomes 


A2(t) 


faj.  (t)  + nsp Ct)  + cos(2irfc|  + *TH(t)}  n15(t) 


+ cos(2tt£c  + <f»3(t)}  nTH (t)  ] c os27rfc(t  - 

+ [sin{2Trfc|+  <t>TH C t^ } n15(t)  + sin{  2trf  c^r-  + $3 (t)  }n^H(t)  ] ■ 


sin2iTf  (t  - — ) 
cv  c ’ 


(170) 


A phase  diagram  of  Eq.  (170)  is  plotted  in  Fig.  17. 

In  Fig.  17  the  desired  phase  to  be  measured  in  2irf 


2 2 
c c 


2z 


Due  to  noise,  the  phase  that  is  actually  measured  is  27rf  — + <J>  , 

CO  v 

where  4>e  is  the  phase  error  or  deviation  due  to  noise. 

The  phase  deviation  is  written  as 


$ = tan 


-1 


n'H(t)Sinei  + n5(t)sin02 
o4H+nsptt)*nTH(t)cos9l*n15(t)cos02 


] (171) 


where 


6i  ■ 2*-r  ¥ + +3ct) 


(172) 


and 


9„  = 2tt 


T c + 4,THCt) 


(173) 


Eq.  (171)  is  not  very  enlightening.  It  becomes  somewhat  clearer 
if  looked  at  in  it  two  limiting  cases:  (1)  speckle  noise  dominant 
and  (2)  quantum  noise  dominant.  Each  of  these  cases  is  discussed 
below. 
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sin2nf 


If  speckle  noise  is  dominant  over  quantum  noise,  then  all 
quantum  noise  terms  can  be  neglected  and  the  phase  deviation  is 
due  to  speckle  noise. 

In  Eq.  (171) , the  n'H(t)  and  n1&(t)  terms  are  due,  at 
least  in  part,  to  quantum  noise.  If  they  are  neglected,  the 
phase  deviation  due  to  speckle  noise  is 


d>  = tan 
sp 


-1 


"sp(t) 


a,2(t) 

xh 


-]  = 


(174) 


From  the  last  equation  it  appears  that  the  phase  of  the  laser 
field  can  be  measured  exactly  for  speckle  noise  limited  systems. 
This  seems  intuitive  because  the  speckle  noise,  like  the  laser 
field,  is  propagated  from  the  surface  to  the  imaging  system. 

Since  both  take  the  same  path,  they  should  suffer  the  same  delay. 
Most  imaging  systems  are  speckle  noise  limited  therefore,  it 
appears  they  should  work  well  for  ranging. 

On  the  other  hand,  if  quantum  noise  is  dominant,  the  n (t) 
and  n^(t)  terms  are  neglected  (n^(t)  results  from  the  cross 
correlation  of  the  speckle  and  quantum  noise,  thus  it  is 
neglected  for  both  cases) . The  phase  deviation  due  to  quantum 
noise  is 


-1 


n'H(t)sin0 


<J>TH  = tan  [ 

iri  a.  (t)+n'Hcos0. 


LTH 


(175) 


This  result  is  very  similar  to  the  result  obtained  by  Ziemer 

and  Tranter  (Ref  24:277-278)  for  the  case  of  angle  demodulation 

in  the  presence  of  noise.  (Eq.  175)  does  not  indicate  the 

magnitude  of  the  phase  error,  but  it  does  imply  that  some 

deviation  occurs.  It  can  be  seen  that  for  a.  (t)  > > n~„  , 

1TH  1 

the  phase  deviation  is  small.  If  the  mean  and  covariance  of 
Eq.  (175)  are  found,  then  the  deviation  can  be  better  understood. 
However,  that  is  beyond  the  scope  of  this. thesis. 

One  practical  consideration  which  should  be  taken  into 
consideration  for  both  cases  is  that  phase  detectors  are  amplitude 
dependent.  All  the  components  which  make  up  the  phase  deviation 
term  (Ref  Eq.  (171))  have  random  amplitudes.  Therefore,  the  results 
of  both  cases  are  true  only  for  a carefully  designed  phase  detector 
which  is  insensitive  to  amplitude  variations  of  the  signal. 

The  important  result  of  this  chapter  is  that  it  appears  that 
laser  line-scan  imaging  systems  can  be  effectively  used  for 
ranging.  The  only  assumptions  used  in  showing  this  were: 

(1)  the  bandwidth  of  the  laser  source  is  narrow  enough  that  laser 
field,  at  any  point  in  space,  is  constant  across  its  width  and 

(2)  the  imaging  system  is  speckle  noise  limited.  Neither 

one  of  these  are  severe  restrictions.  Also  important  is  Fig.  13. 

It  gives  the  values  of  fTC  and  f that  can  be  used  to  avoid  aliasing 
This  insures  that  the  correct  image  is  processed  by  the  system. 
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VI . Conclusion 

Summa  ->v 
— 

TS  cts  of  speckle  on  the  output  current  of  an  imaging 

sy  e examined.  The  imaging  system  was  composed  on  a laser 

J 

scuu  xi. 0 a rough  surface  (compared  to  the  laser  wavelength) 
and  a receiver  measuring  the  field  reflected  from  the  surface. 

Three  types  of  lasers  and  two  types  of  receivers  were  examined. 

The  lasers  were  designated  by  their  bandwidths  as  monochromatic, 
quasimonochromatic , and  broadband.  The  two  receivers  were  the 
direct  detection  and  heterodyne. 

The  field  reflected  from  the  rough  surface  was  "crudely" 
modeled  by  multiplying  the  incident  laser  field  by  a reflectance 
term  and  a random  phase  term.  The  reflectance  term  was  the 
desired  quantity  or  "signal"  to  be  measured.  The  phase  term 
simulated  the  effects  of  speckle.  All  fields  were  propagated 
by  the  Huygens -Fresnel  integral. 

The  system  composed  of  a monochromatic  laser  and  a direct 
detection  receiver  was  examined  first.  A system  function  dependent 
on  the  laser  beam  spot  size  at  the  surface  was  defined.  It 
was  shown  that  the  resolution  ability  of  the  system  was  determined 
by  the  width  of  the  system  function.  By  second  moment  techniques, 
the  "noise"  (due  to  speckle  effects)  measured  by  the  system 
was  found  to  be  signal  dependent.  The  significance  of  this 
result  was  that  the  signal  to  noise  ratio  (SNR)  could  not  be 
improved  by  increasing  the  signal  power.  The  value  of  the  SNR 
was  shown  to  be  greater  than ‘or  equal  to  one.  The  far  field 
and  near  field  cases  were  found  to  be  identical.  The 
monochromatic  laser  was  then  replaced  by  a quasimonochromatic  one. 


The  results  for  this  new  system  were  shown  to  he  identical  to 
the  monochromatic  results. 

For  the  imaging  system  consisting  of  a heterodyne  receiver 
and  a monochromatic  laser,  the  measured "noise"  due  to  the  effects 
of  speckle  was  also  found  to  be  signal  dependent.  Thus  the 
SNR,  which  was  equal  to  one,  could  not  be  improved  by  increasing 
the  signal  power.  The  resolution  of  the  system  was  shown  to 
vary  with  the  beam  spot  size  at  the  surface  and  the  spatial 
fourier  transform  of  the  receiver's  detection  area.  The  far 
and  near  field  cases  were  identical.  The  heterodyne  results 
were  unchanged  when  the  monochromatic  laser  was  replaced  by  a 
quasimonochromatic  one.  When  a broadband  laser  was  used,  the 
only  result  which  changed  was  that  the  resolution  ability 
of  the  heterodyne  system  improved. 

The  most  significant  results  occurred  when  the  imaging 
system  was  composed  of  a direct  detection  receiver  and  a broadband 
laser.  It  was  shown  that  this  system  is  not  effected  by 
speckle  noise.  This  was  different  from  all  other  systems. 

They  were  limited  by  speckle  effects.  This  system  was  found  to 
have  the  best  resolution  characteristics  of  all  the  systems 
investigated.  From  these  results,  it  was  determined  t^at  this 
system  is  the  optimum  system  to  use  for  laser  line-scan  imaging. 

As  an  application,  the  system  consisting  of  a quasimonochromatic 
laser  and  heterodyne  receiver  was  used  for  ranging.  In  addition 
to  the  effects  of  speckle,  quantum  effects  were  included.  The 
laser  field  was  amplitude  modulated.  The  phase  delay  in  the 
modulation  term  of  the  reflected  field  was  measured  and 


used  to  determine  the  distance  from  the  surface  to  the  system. 
However,  the  measured  phase  delay  included  a phase  error  due  to 
quantum  and  speckle  effects.  It  was  shown  that  for  the  case  of 
speckle  noise  dominate  over  quantum  noise,  the  phase  delay  could 
be  accurately  measured.  For  the  case  of  quantum  noise  dominate, 
it  was  found  that  an  error  would  result  in  the  measured  phase 
delay.  No  attempt  was  made  to  quantify  this  error.  It  was 
noted  that  the  obtained  results  were  only  valid  for  a phase 
detector  which  is  insensitive  to  amplitude  variations.  It 
was  concluded  that  ranging  could  be  effectively  performed.  A 
method  for  selecting  the  intermediate  and  carrier  frequencies 
used  in  the  ranging  problem  was  also  presented. 

Recommendations 

In  this  thesis,  a laser  line-scan  imaging  system  was  modeled 
for  different  combinations  of  laser  sources  and  receivers.  These 
models  were  based  on  theoretical  analysis.  It  is  recommended 
that  these  models  be  compared  to  experimental  data  to  test  their 
validity. 

One  calculation  used  to  determine  the  effects  of  speckle 
on  broadband  imaging  systems  (Ref  Eq.  (123))  was  only  an 
approximation.  It  seems  feasible  that  this  calculation  can  be 
quantified  by  comparison  to  the  analysis  of  similar  calculations. 
This  would  add  to  the  understanding  of  the  broadband  results  in 
this  paper. 

In  the  ranging  problem,  the  phase  error  was  presented  but 
not  analyzed.  A statistical  analysis  of  the  phase  error  would 
result  in  a better  understanding  in  the  magnitude  of  the  error 
and  the  variations  in  its  amplitude.  As  previously  stated, 
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phase  detectors  are  amplitude  dependent.  A statistical  analysis 
would  indicate  the  effect  of  amplitude  variations  on  phase 
detectors.  Also,  the  ranging  problem  was  examined  for  the  case 
of  a constant  surface  reflection  coefficient.  If  the  problem 
were  solved  for  a spatially  varying  coefficient,  a more  general 
solution  would  be  obtained. 
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Appendix  A 


Extended  Huygens-Fresnel  Integral 

Recall  from  the  section  on  complex  representation  that  the 
real  scalar  optical  field  is 

u(r,t)  = ACr.t)  cos  [2irr  t-<{>  (r  ,t)  ] CA-1) 

and  its  complex  envelope  (with  exp  [- j 2TrfQt]  time  dependence) 
is 


U(r , t)  = ACr,t)  exp [ j $(r ,t) ] 


CA-2) 


where  r is  a vector  with  spatial  coordinates  (x,y).  Let 
the  temporal  Fourier  transform  of  Eqs.  (A-l)  and  (A-2)  be 
denoted 


Ft[uOF,t)]  = VCr  ,f) 


(A- 3) 


and 


Ft[U(r,t)]  = VLp(r,f) 


(A-4) 


By  Fourier  analysis  techniques  it  can  be  shown  that  Eq.  (A- 3)  can 
be  written  in  terms  of  Eq.  CA-4)  as 

V(r,f)  = ^Ft[U*(r,t)  exp  (j  2irfQt)  ] 


+ j Ft[U(r,t)  exp(-j27rfot)] 


(A-5) 
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For  the  case  where  the  complex  envelope,  U(r,t),  is  monochromatic, 
the  Fourier  transform  of  the  optic  field  can  be  written 

. V(r,f)  *=  V(r,fQ)«(f-f0)  ♦ V(r,-f0)«(f+f0)  (A-6) 

The  amplitude  of  the  last  term  in  Eq.  (A-6)  is 
vCr,-f0)  - \ Ft[UCr,t)exp(-j2irfot)] 


VrDCr,f+fo) 


I VLP 


(A-  7) 


The  two  dimensional  form  of  the  Huygens-Fresnel  integral 


is 


U fr  1 

ulLrlJ  jt 


exp [ j kz ] 


U (r) exp [ j lr1'rl2]dr 


C A - 8 ) 


where  X=  / £ . On  comparing  Eqs.  (A-5)  and  (A-7)  it  is  seen 
o 

that  Eq.  (A-8)  can  be  written  as 


Vi  Of  !•-«.)  - s^1 


V(r,-fo)exp[j^-|r1-r  j 2]dr  (A-9) 


After  completing  a change  of  variables  f*-f  , Eq. 


Vl(rl»f)  = 


exp[-j2rrz- 

-jfz 


V(r,f)cxp[-j^-|r1 


(A-9)  becomes 
r 1 2 ] d v (A- 10) 


This  equation  is  valid  for  propagation  of  any  field  component 
at  frequency  f.  To  express  the  results  in  terms  of  complex 
envelopes,  note  that  for  all  frequencies  f near  -fQ 


VCr.f)  - £vLp(7,f  ♦ ffl) 


(A- 11) 
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Substituting  Eq.  (A-ll)  into  Eq.  (A-10)  yields 


Vl,LP(rl,f  + fo3 


£ exp  [ - j 2ttz—  ] 


JCZ 


exp[-j^|  |r1-r|2]dr 


(A- 12) 


The  temporal  inverse  fourier  transform  of  Eq.  (A- 12)  is 

2 


U^r.Oe 


^nf  t 


2nzc 


k , f HT  {u(r»t"l  ' ' 2zc'  ) 


r,  -r 


exp[-j2irfo(t-|-  — )])dr  (A-13) 


By  use  of  the  chain  rule,  Eq.  (A-13)  becomes 


ui'ri-t>  ■ life  exp[jkl1 


exp  [^fei  r1‘ri  *■  ] 


-'-fc-  )-j2»f0UCr,t-|-  -4^  ))dr  (A-14) 


Ir.-rl2  _ |r,-r|2 


Eq.  (A-14)  can  be  simplified  somewhat.  If  it  is  assumed 
that  the  bandwidth  of  U(r,t)  is  much  less  than  the  optical 
frequency  fQ,  then  the  f in  the  numeration  of  Eq.  (A-12) 
can  be  approximated  by  fQ.  The  resulting  temporal  inverse 


Fourier  transform  of  Eq.  (A-12)  is  then 


. m2 


Vr>t)  * ^xk-1  JJ  Ul(r»t‘l"  "4^  )cxp[j2\|r1-r|2]dr 


(A- IS) 
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This  assumption  implies  that  the  field  does  not  change 
significantly  in  time  relative  to  the  variation  at  fQ. 

Eq.  CA-15)  would  also  result  is  the  derivative  term  in  Eq.  (A-14) 
was  assumed  to  be  insignificant  in  comparison  to  the  other  term 
in  that  equation. 

Another  simplification  can  be  made  to  Eq.  (A-14)  when  the 
bandwidth  of  U(r,t)  is  restricted  by 


1 _ _ ^01 1 max 
F > > c 


(A-16) 


where  is  defined  by  Eq.  (10).  This  is  the  quasimonochromatic 
restriction  (Ref  12:108).  The  resulting  Huygens-Fresnel 
integral,  which  is  valid  for  quasimonochromatic  fields, 
is 


UiCFj.t)  = C-X^-k-Zl  J j U(r,t)exp [j jy  |7rr|2]dr  (A-17) 

The  Huygens-Fresnel  integral,  which  is  valid  for  monochromatic 
fields,  has  been  extended  so  that  it  is  valid  for  fields  of 
any  bandwidth  (Ref  Eq.  (A-14)).  Two  special  cases,  the  broad- 
band and  quasimonochromatic,  were  presented.  The  propagation 
of  broadband  fields  is  given  by  Eq.  (A-15)  and  the  propagation  of 
quasimonochromatic  fields  is  given  by  Eq.  (A-17). 
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Appendix  B 

The  Method  of  Stationary  Phase 


The  method  of  stationary  phase  is  an  asymptotical  approximation 
to  the  solution  of  an  integral.  It  is  discussed  by  Papoulis 
in  two  of  his  books  (Ref  13:  234-250,  15:139-143).  The  method 
is  summarized  below. 

Consider  the  integral 

Q(Af)  - | | gx  (ot , B) exp  [ j Af  g2(a,B)]dadB  (B-l) 

R 

where  R is  a region  in  the  a,B  plane,  g 2 C01  * & D must  be 
twice  differentiable  in  R.  The  extreme  value  of  g2(ct,B) 
at  a single  point  Ca0»&0)  in  R is  found  by  the  first  partial 
derivatives  of  g2 C°t , B ) : 


3g2(a,B) 

da 

3g2(a,8) 
tb — 


■ 0 ■ fW'W 

0 


If  g j (a , 8)  is  continuous  at  (aQ,Bo)  and  if 
g2,aag2,BB  '(g2,a6)  * 0 

and 


(B-  2) 


(B-3) 


(B-4) 


g2,  BB 


t .0 


CB-5) 
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where  g2 
evaluated 

Q(flf) 


, g,  , and  g,  are  the  partial  derivatives 
»aa  * * »ag 

at  Cao»0o)»  then  as  Af  approaches  infinity 


j2vg1(ao>go) 

^g2,aag2,6g"  g2ta^^ 


exp[jAfg2Ca0,60)] 


(B-6) 
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